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ABSTRACT 


Prismatic and concentric helical prismatic dislocations are generated by the 
stress field which appears around a spherical glass inclusion in a silver chloride 
erystal when the crystal is cooled from a temperature of 370°C to room tem- 
perature. The properties of these dislocations are described and their origin 
is discussed. 


§ 1. INTRODUCTION 


DuRrtincG the course of experimental work on chemical sensitization (Evans 
et al. 1955), it was noticed that the chemical and photochemical reactivity 
of silver halide crystals was enhanced in the neighbourhood of glass 
particles. Chemical development was often spontaneously initiated near 
the particles and this source of fog was eliminated when they were removed 
by filtering the molten silver halide through fine glass capillary tubes. 
At the time, it was suspected that the particles increased the average 
local density of dislocations and that this was responsible for the increased 
reactivity. We have now deliberately incorporated particles of hysil 
glass in silver chloride crystals and studied the distribution of dislocations 
around them by etching (Jones and Mitchell 1957) and by the separation 
of silver along the dislocation lines during exposure (Hedges and Mitchell 
1953). Prismatic (Seitz 1950) and helical prismatic dislocations with 
{110) directions as axes are generated by the stress field which arises on 
cooling from the differential contraction between the silver chloride and 


+ Communicated by the Authors. 


P.M. A 


2 D. A. Jones and J. W. Mitchell on 


“a 


the glass. Spherical particles have also been made and introduced into 
the crystals ; the symmetry of the resulting stress field has greatly facili- 
tated the interpretation of the observations. 


§ 2, EXPERIMENTAL METHODS 


The silver chloride was prepared and purified as described by Clark and 
Mitchell (1956). Glass spheres were made by the method of Sollner (1939). 
Glass particles were introduced into the oxygen stream of an oxy-coal gas 

_blowpipe and passed through the jet where they fused into spheres. These 
were collected in a dish containing distilled water and separated into size 
classes by sedimentation. The most useful particles were of almost perfect 
spherical form with diameters between 0-5 and 5 microns. 

Single crystals of silver chloride were made in the form of thin sheets by 
the method described by Clark and Mitchell (1956). A drop of a suspension 
of the glass spheres in distilled water was placed on the surface of the upper 
glass plate which was to be in contact with the silver halide. After the 
evaporation of the water, a reasonably uniform distribution of the spheres 
remained. Their surface density could be controlled conveniently by 
varying the concentration of the suspension so that isolated or interacting 
prismatic dislocations could be produced as desired. The spheres became 
imbedded in the crystals but usually remained within 20 microns of the 
surfaces ; this was convenient for etching, for the separation of silver along 
the dislocation lines during exposure and for microscopic examination. 
After cooling to room temperature, the crystals were separated from the 
glass plates, between which they were grown, and cut into sections. These 
were annealed in an atmosphere of chlorine for 8 hours at 370°c. The 
chlorine was then condensed by immersing a side arm in liquid air and the 
specimens were cooled to room temperature. This treatment reduces the 
density of dislocations in many crystals with (001) surfaces to a negligible 
value. The properties of the prismatic dislocations which are generated 
around spherical inclusions during the cooling of the annealed crystals can 
therefore be studied without complications arising from their interactions 
with other dislocations. 

In the early stages of the work, the crystals were etched and then 
exposed to render visible the dislocation lines below the surface. It was 
found, however, that etching, which removed the upper half of the loop 
systems by dissolving away the surface, destroyed information which was 
essential for the interpretation of the observations. The crystals have 
since been simply exposed to a 250 watt medium pressure mercury vapour 
lamp with a frosted envelope at a distance of 10 em for 2 to 3 minutes 
During this time, very little silver separates on (001) surfaces to obaeut 
the dislocations which are decorated down to a depth of 20 to 30 mriceee 
Tiong thedlilovstion ere aatoeentt ie eee ee 
can searcely be resolved Sane ae a pene DAS 
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decorated than edge dislocations. The paper is illustrated with photo- 
graphs of wire models of the dislocation systems. The illustrations are 
effectively projections of the systems on an (001) plane from a point above 
the plane. 


§ 3. EXPERIMENTAL OBSERVATIONS 


Dislocation loops are formed around any glass inclusion in a crystal of 
silver chloride which has been grown from the fused substance. The 
observations which are described in this section have, however, all been 
made with spherical particles. 

The thin single crystal sheets of silver chloride usually have (001) surfaces. 
The systems of prismatic and helical prismatic dislocations which are 
formed around a sphere have as possible axes the four (110) directions 
which lie in a plane passing through the centre of the sphere parallel to 
the surface of the crystal. One, two, three and four systems have been 
observed around individual spheres. Sometimes a few loops of systems 
along the ¢110) axes which are inclined downwards from the surface have 
been observed. In crystals with (111) surfaces, systems of prismatic 
dislocations are formed with any of the six possible (110) directions. 
which are parallel to the surface as axes. The systems formed at different: 
particles along particular (110) directions in a plane parallel to the surface 
will be illustrated and discussed in the remainder of this section. 

The majority of the systems consist of sequences of loops and figures of 
eight as shown in figs. 4 and 7, together with more complicated closed 
dislocation loops such as those illustrated in figs. 5 and 6. The diameter 
of the sphere in fig. 4 is 3-5 microns. The diameter of the loops is 2-4 
microns and the overall length of the system 15 microns. The loops have 
a nearly constant diameter over the whole length of a system and this is 
usually slightly smaller than the diameter of the spherical inclusion. 
When they are far from other inclusions and systems of prismatic dis- 
locations, they usually lie in planes normal to the axes of the systems along 
which they are spaced at remarkably regular intervals. 

In figures of eight and more complicated configurations (fig. 6), the 
dislocations cross at an angle of approximately 90°. Although the 
diameters of the spherical particles were usually about 5 microns and the 
diameters of the loop systems slightly less than this, the quality of the 
decoration allows the details of the configurations to be clearly resolved. 
It has been found that the complex configurations are formed from two 
concentric helical dislocations of opposite twist. At all the cross-over 
points, one helical dislocation is outside the other, showing that they must 
have been generated as concentric helices, and not as two independent 
helices which were subsequently combined by lateral displacement. 

In the many hundreds of systems which were examined in detail, it 
proved possible to find configurations which appeared to be in every 
possible stage of development. Figures 1 and 2 show configurations 
which were frequently observed. It was deduced that a narrow loop was 
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first formed in the silver chloride at its interface with the sphere by glide 
in a plane containing a (110) direction and that it expanded around a 
cylindrical surface with the (110) direction as axis to generate single 
loops, figures of eight and more complex configurations. These units were 
randomly distributed along the (110) directions but a concentric helical 
dislocation was often observed in contact with the sphere as in fig. 3. 
Here the diameter of the spherical particle was 4-1 microns and the 
diameter of the prismatic dislocation 4:1 microns near the surface of the 
sphere, decreasing to 3-3 microns at the other end of the system. 
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Many single helical dislocations with axes along (110) directions were 
observed. The diameter and pitch usually varied with distance from the 
sphere and the dislocations could often be followed to the point where 
they made contact with the surface of the crystal. Single prismatic 
loops around single helical dislocations, around the inner helical dis- 
location of a concentric pair, and around both dislocations of a concentric 
pair have been observed. Such isolated prismatic dislocations have never 
been observed to surround only the outer dislocation of a concentric pair 

A number of specimens were prepared with a high density of Aare 
so that the interactions between the prismatic dislocations could 7 
uid In these specimens, loop systems generated by one sphere were 
sometimes observed to. pass by a smaller particle which had not generated 
any dislocations. The individual loops were then displaced as shown in 
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Fig. 
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Fig. 
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fig. 7. A similar displacement was produced if a loop system passed by 
another loop system with an axis at right angles and it became evident 
that the inclination of the prismatic dislocations provided a sensitive: 
indicator of the presence of local shearing stresses in the crystals. 


§ 4. INTERPRETATION OF THE OBSERVATIONS 


The interpretation of the observations is based upon an extension of a. 
mechanism proposed by Seitz (1950) in his discussion of prismatic punching 
and prismatic dislocations in thallous halide crystals. He pointed out 
that two prismatic dislocations of opposite sign could be generated by the 
expansion of a dislocation loop around the glide prism. Nye (1949) has. 
established that silver chloride deforms by pencil glide with a (110) slip 
vector so that prismatic punching is a possible mode of deformation. 

To interpret our own observations, we have to explain the generation 
of prismatic and helical prismatic dislocations on the surfaces of cylinders. 
with axes along the possible ¢110) slip vectors and diameters slightly 
less than that of the spherical inclusion. The motion of such prismatic 
dislocations along the slip cylinder transports material away from the 
spherical inclusion to reduce the strain arising from the differential con- 
traction on cooling. It may reasonably be assumed that there is no 
coherence between the silver halide and the glass at the interface and that 
the silver chloride there is redistributed by subsidiary slip and diffusion. 
processes. The prismatic dislocations move outwards without any appre- 
ciable change in their projected area upon a plane at right angles to the 
axis of the glide cylinder; diffusion processes leading to climb are not 
important for these phenomena. The experimental observations show 
that the strains are transmitted over distances large compared with the 
diameters of the glass spheres. 

The stress field produced by differential contraction between the silver 
chloride and the glass sphere has, ideally, spherical symmetry. The 
shearing stress is zero in radial directions. With respect to any (110) slip. 
vector, it has a maximum value at the surface of the sphere upon a cylin- 
drical surface which has the slip vector as axis and a diameter 1/2 times. 
the radius of the sphere. Small dislocation loops are therefore formed on 
the surface of this glide cylinder. The forward edge component then 
glides outward under the influence of the steadily decreasing stress field. 
The screw components of the loop, which are parallel to the axis of the 
glide cylinder, experience a tangential force which causes them to glide 
in opposite directions around its surface. This force is greater near the 
interface than further away so that, although the dislocation must retain 
the screw orientation at the interface, it will elsewhere become inclined to 
the axis as it glides round the cylinder. The introduction of the edge 
components as a result of this inclination stabilizes the dislocation on a 
particular glide cylinder because any radial displacement would involve 
climb. If these two dislocations remain upon the same cylindrical 
surface during their rotation round the axis under the action of the 
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shearing stress, their meeting will produce a prismatic dislocation corre- 
sponding to the introduction of additional material into the lattice. It 
is convenient to describe this as a positive prismatic dislocation. An 
equivalent amount of material is removed from the incoherent interface 
between the silver chloride and the glass sphere. A succession of such 
dislocations may be generated by independent events in which small 
loops are formed by slip and expand round glide cylinders. The 
remarkably uniform diameter of the loops in many systems is probably 
determined by the diameter of the glide cylinder upon which the shearing 
stress at the interface has the maximum value. The regular spacing 
arises from their mutual repulsion which is also responsible for the trans- 
mission of stresses along the cylinder. This is essentially one of the 
mechanisms proposed by Seitz (1950) for the generation of prismatic 
dislocations and for the interpretation of prismatic punching. 

The observations illustrated in fig. 7 confirm that positive prismatic 
dislocations are formed. The prismatic dislocations pass another spherical 
inclusion which is surrounded by a stress field. The strain energy is 
lowered by the overlapping of the region of dilatation around the edge of 
the positive prismatic dislocation and the region of compression around 
the spherical inclusion. This interaction, together with an analysis of the 
forces acting on the dislocations, establishes that they are of a positive 
type. 

The most striking feature of the observed phenomena, the occurrence of 
figures of eight, of more complex closed dislocation loops (fig. 6) and of 
concentric helical dislocations (fig. 3), requires an extension of the above 
mechanism for its interpretation. If the dislocations are displaced 
radially during their rotation round the glide cylinder, they may not meet 
to generate prismatic dislocations, but may continue their rotation with 
the production of concentric helical dislocations. This is illustrated by 
the loops in contact with the spheres in figs. 1, 2 and 3. If the dis- 
locations meet after each has made a complete rotation around the glide 
cylinder, figures of eight will be formed. The more complex systems 
result from a larger number of rotations before intersection. The radial 
displacements which lead to the formation of helical prismatic dis- 
locations and to their subsequent intersection are most likely to occur 
where the dislocations involved are in the screw orientation, that is to say, 
at the interface with the spherical inclusion. As already pointed out, 
when climb is not permitted, the edge components of composite dis- 
locations stabilize them against radial displacements with respect to the 
axis of the glide cylinder. This also explains why the dislocations do not 
attract each other when they are rotating on their respective clide 
cylinders in the same way as two screw dislocations of opposite hand 
The dislocation lines cross at an angle of nearly 90°. The attractions of 
aaa pectneats cae then, in a sense, be balanced by the mutual 

s g ponents although a displacement of the edge 
components could only result from climb. : 
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In the majority of the concentric helical prismatic dislocations, the 
outer and the inner helices were of very nearly the same diameter but in 
some (as, for example, in fig. 3), the outer helix was of constant diameter 
while the diameter of the inner helix decreased as the surface of the sphere 
was approached. In these cases, the screw element of the inner dis- 
location at its intersection with the interface has moved inwards during 
its rotation round the axis of the glide cylinder. 

The single helical dislocations are formed when the primary loop 
expands so that one of the screw dislocations meets the surface of the 
crystal. The rotation of the screw element of the other dislocation at its 
contact with the interface between the silver halide and the spherical 
particle around the axis of the glide cylinder will generate a single helical 
dislocation. 

The single prismatic dislocation loops which surround the helical dis- 
locations are produced by the expansion of loops, formed by slip, around 
cylindrical surfaces with greater or smaller radii than the outer dis- 
location of the helical system. 


§ 5. CONCLUSION 

A stress operated mechanism for the generation of helical prismatic 
dislocations has not previously been described. The method of in- 
corporating glass spheres into silver halide crystals, which has been intro- 
duced for the first time in this investigation, provides a powerful means 
for generating helical prismatic dislocations and for studying their 
properties. 
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The Centre of a Dislocation: II—The Dilated Slit} 


By EvizaBetu H. YOFFE 
Cavendish Laboratory, Cambridge 


[Received October 14, 1957] 


ABSTRACT 


This paper describes the second part of some investigations of the elastic 
field near the centre of an edge dislocation. It is concerned with the local 
expansion due to atomic misfit, causing a singularity similar to a particular 
type of centre of pressure. The corresponding stress field is added to that 
of Part I to give the complete elastic field of a dislocation in a crystal. 

The solution is also applicable to the dilation of a slit by a solid inclusion, 
and therefore gives an indication of the stress field of a tapered precipitate 
in a crystal matrix. 


§ 1. [INTRODUCTION 


Av the centre of an edge dislocation in a crystal lattice the atoms next 
the slip plane are so far out of position that they are necessarily displaced 
normal to the slip plane as well as parallel to it. The atoms as it were 
ride up over those on the opposite side of the slip plane in the limited 
region where the relative transverse displacement is large (cf. Read 1953). 
This has the effect of introducing a discontinuity in the elastic field, as 
though additional elastic material had been inserted in a finite slit along 
the slip plane. 


Fig. 1 
yp 
re 
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The dilated slit and Cartesian axes. 


The problem considered in this paper is that of determining the stress 
field due to such dilation of a slit, for various distributions of the eddicenn 
material. Similar boundary value problems have been solved b 
Muskhelisvili (1953) using the conformal transformation of an ellipse i 
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a circle. The slit is treated as the limit of a narrow ellipse and a solution 
is given for the straight slit dilated by uniform internal pressure to an 
elliptical aperture. 

Since the final form of the slit in the present problem is not necessarily 
elliptical, it is found more convenient to use Cartesian coordinates, and a 
different method based on higher order dislocations, as described in a 
recent paper by the author (Yoffe 1957). 


§ 2. Dination Dus To Miserrr 


For purposes of calculation the metal crystal containing a dislocation 
is replaced by an isotropic elastic body in a state of plane strain, and 
Jartesian axes are so chosen that the z axis represents the dislocation 
line and the (vz) plane corresponds to the slip plane. The problem is to 
determine the components of displacement u,v, in the a, y directions 
respectively, and the stress components 22, yy, xy due to various dilations 
of a slit between the points (h, 0) and (—A, 0). 


Hiey2, 


The quantity d in a two dimensional lattice. 


The amount of dilation, or relative displacement wv’ of the two sides 
of the slit, is defined as a function of x on the slip plane : 
0 (x)= lim v(x, €)—v(a, —e) 
e>0) 
and is assumed zero for | w |>/, increasing symmetrically to a maximum 
value d at z=0. This is because v’ is due to atomic misfit caused by an 
edge dislocation of relative displacement wu’, where w’ increases from zero 
at «=h to one atomic spacing a at x=—h. 
In the simple two dimensional array sketched in fig. 2 the maximum 
value d is given by : 
d=a(1—cos 30°)=0-134a. 
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In more general three dimensional lattices 
a—Co 
where c is a constant to be calculated in each case and is of the order 0-1. 


If the atoms are regarded as spheres rolling up over one another as a 
fig. 2, the following simple relation between w’ and v’ may be obtained : 


d : : fs 4 ay —cos7 | 
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but w’ must still be determined as a function of x by choosing a particular 
Law of Force (Yoffe loc. cit.). 


Fig. 3 


aaa 2h SSS 


Alternative forms assumed for w’(z). 


It is assumed that v’(a) resembles one of the following curves, represented 
for simplicity in terms of straight lines and parabolas (see fig. 3) : 

(a) Straight line of arbitrary initial gradient merging into a parabola 
(A— Bz?) near x=0. 

(6) Parabolas at x=0, x=-¢h, joined by straight line segments. 


When v’ is assumed to take one of these forms the stress field is readily 
derived from second order dislocations. 

The final displacements of the atom centres are not assumed but 
calculated. The quantity v’ represents the width of imaginary additional 
material inserted in the slit. It may be considered as the non-elastic 
part of the relative displacement of the atoms nearest the slip plane. 
The final displacement Vv contains both elastic and non-elastic terms. 


§ 3. DisLocaTions OF THE First THREE ORDERS 


All the dislocations used in this paper are plane dislocations of the 
second type, requiring imaginary additional material along the negative 
«axis. The order of a dislocation depends on the variation of v’ with x. 

The zero order dislocation of this type has v’=constant on «<0. It 
has the stress function 


Uj=Aalogr... (2 =a +4?) Lethe eS) 
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for some constant A, from which may be deduced the components of 
stress and displacement : 


Wo=(Aax|r*)(3—(2x%/r?)) i 
2pvy=A[A(1—o)0—(ey/r)] 
where tan@=y/x. This is a Volterra dislocation, and corresponds to a 
Burgers edge dislocation on the (yz) plane. 

Vp =27(1—o)A/p. 


The first order dislocation is also a Volterra dislocation, since it has 
relative displacement as a rigid body : 


yoccronz<0. 


A stress function having the required properties is 


U,=—}A(2r? log r—r?) eee hae (8) 
with components of stress and displacement such as : 
YY = —A (log r+ cos? 6), ] 
2uv,=—A[(1—2c)y log r+ 2(1—o)20— 2(1—o)y], (4) 
V4'=2n(1—oc)Ar|p. 


‘This dislocation corresponds to a linear distribution of zero order disloca- 

tions along the negative x axis (cf. Mann 1949). The zero order stress 

function may therefore be obtained by differentiating U, with respect 
to —x: 2 

O x wi eu ee 
0x 


.and the same is true of the various components 


Uy 


— a =H¥p ete. 
(cf. Yoffe 1957). e 
The second order dislocation is not a Volterra dislocation since it haa 
relative displacement proportional to the square of the distance from the 
origin. A suitable stress function is 
1 2x3 + 5 7 
(ome? UMN 5.52 op eh ah? ee ae) 
v= zal ( 5 +2ny*) log 3 V0 9% zy | (5) 
with components : 
£2,—=A (ax log r—2y0—x), 
YY,—=Ax log r, 
Ly y= — Ay log r, 
2uu =A {2 log r[(1—2c)x?—(3—2c)y?|—8(1—o)ay8 —¢ (6) 
—(3—40)x?-+(5—4o)y?}, 
2uv,= A[(1—20)ay log r+ (1—o)2?0-+ oy?8— (1— 20) xy], 
Vv, =7(1—o) Aru. 
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It is readily verified that 
OU 5 es ee ee 


with similar relations for each set of components. 

It may be noted that all terms of (6) are finite at the origin, so that a. 
continuous stress field may be built up by combining second order disloca- 
tions. Such a field would however have singularities at infinity unless. 
the dislocations are so chosen that their fields cancel at large distances. 

By an argument developed in the previous paper (Yoffe loc. cit.) this 
cancelling is achieved by finite difference approximation to the derivative. 
Let each term of eqn. (7) be differentiated once more with respect to 


NEE aU, OU,  ay_,, : 
Giese ues 3. Che ee (8) 
then the stress function U_, gives a suitable stress field at infinity : 
U_,=—A (log r+ cos? 9), ole, tee i RO 
ae A Say? 
= 5(1— a te | 
A Lae ee | 
C= 5 (3 4 =) ; 
= Ax Ba4 
LY = a (o =) ; Me (10), 
. AZ x2 — 3y? 
2uu_ y= sa | 4 )+ 2 ih 
ip ee 
2,0 == a 5—4o+ y i ] 


es h ces a 
U_, is the stress function of a centre of pressure at the origin, causing: 
compression in the y direction and tension in the « direction. All. 
components of U_, tend to zero at infinity. 


Fig. 4 


ye =) € az " t ph 7 
The relative displacement v’(x) obtained from two zero order dislocations 


. If therefore dislocations of zero, first or second order are combined’ 
mn such a way as to represent one of the terms of eqn. ( 8), then there will 
be no singularity at infinity. | 
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For example it is readily verified that a zero order dislocation at 
(h, 0) and an equal but negative one at (—h, 0) may cause the slit to be 
dilated by a constant amount as in fig. 4. Let suffixes h, —h relate 
coordinates to these points as new origins so that 


6,=tan-"(y/x—h), 


r= Vile+h)+y2] ete. 
Then the stress field of the two dislocations may be calculated from the 


stress function U=A[(a—h) log r,—(x+h) log r,] Beis eh (L TP) 


with A=pd/27(1—c). At points a large but finite distance from the 
origin this function U may be expanded in a rapidly converging Taylor 


series : U=U,(«—h)—U,(ax+h) 
ri OU, hF AU, 
Sate ities eta Ox? Seg 


so that U/2h—> —dU,/dx=U_, as in (8), or 
5 dh 2 ibe 
U—> a ar (IGG =CO8* 0) 0 me sa) ete 3 (LS) 
and the field is regular at infinity. 


Fig. 5 


(b) 


(a) Relative displacement corresponding to eqn. (13). 
(b) Relative displacement corresponding to eqn. (14). 
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Similarly the stress fields of first order dislocations may be combined 
in any manner which approximates to their second derivative. For 
example eu, 


U(e—h)—2U (w) + Uy(a-+h) > BS =WU_,. (13) 


ae 
So all components are finite at large distances if first order dislocations 
are placed at (h, 0) (—h, 0) and two negative ones at 0. The resulting 
dilation of the slit is shown in fig. 5 (a), when the constant A in each 
U, is chosen as A =pd/2a(1—o)h. 

A more general form is given by the stress function 


U=U ,(a—h)—U ,(a—k)—U ,(a@+hk)4+ U(x +h) > (h?—k?)U_, (14) 


as shown in fig. 5 (b), where 0<k<h. Equation (14) is readily verified 
by expansion in Taylor series. The constant A in each U, is here 
ued/2m(1—o)(h—k) so that at large distances 


— pa(h+k) 


Cera (oan 


. (log r+cos? 6). Pe pe EY) 


§ 4. Stress FIELD OF THE DiiaTED SLIT 


To dilate a slit along the smooth curves shown in fig. 3, second order 
dislocations are distributed between +/ to give systems equivalent to 
the first term of eqn. (8). They may be distributed in any manner 
representing an approximation to the third derivative, so that at large 
distances the stress function still reduces to CU_, where C is a constant. 

A curve of the type (b) of fig. 3 may be specified by two parameters 
k and k’ which determine the radii of curvature at 0, 2 and —h as follows : 


(Omer (cds) 


d (h—x)? write 
Qh (h—k<av<h) 
y= 4 d (h—x)*— (h—k—2a)? a 
(2h—k—k’) <a ee (k’ <a <h—k) (16) 
eek (h—ax)®—(h—k—a)? _ (k'—a)? ect: 
| (2h—k—k’) ip Shae aoe Cla 


with v’(z)=v'(—a). The dislocations are all of second order, situated 


at the points x=+h, +(h—k), +k’. The stress function of the whole 
system is given by: 


py Valt—h)—U eh +k)  Us(e—k')—U(a+k’) 
k k’ 
U(x«+h—k U Lf, 
je A : ‘ (17) 
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with U, defined by eqn. (5) in which A=pd/n(1—o)(2h—k—k’') and at 
large distances 
U — 3(3h?—3hk-+k®—k2)U_, 


pd (3h2—3hk--k2—k'2) 


ei RES a oy 
=CU_, 


as required. 

A curve of the type (a) of fig. 3 may be regarded as the limit of type (6) 
as k—> 0, the first and last terms of U in (17) reducing to the stress 
functions of first order dislocations : 


U(x—k’)—U (ak! Sh2—k’? 
ee Se Uy(e-+h) >= U_,. 


(18) 


This relation may also be obtained directly, since the curve (@) may be 
specified by one parameter k’ as follows : 


fo (|x| >m 
yi [aia )R—ER) (<r <h) a 
| dh—ae— (li —a)/2K\/(h—4k') (—k'<a<k’) 


| d(h—a)/(hR—4k’) (—hxa<—k’). 


From consideration of eqns. (4) and (6) it is found. that such a v’ may 
be composed of first order dislocations at --4 and second order ones at 
+k’, so the result (18) follows, with A=pd/27(1—o)(h—Hk’). 

If k’ also tends to zero in eqn. (18) the stress system is that of fig. 5 (a) 
and eqn. (14). 

The stress field is completely determined from the stress function, so 
that the problem is now solved. From any given curve of the types in 
fig. 3 the constants k and k’ are determined, substituted in (17) and the 
required components calculated in the usual way. 

For example, the normal stress 7 at points along the slip plane is 
calculated from 


In the case of equally spaced dislocations, k=2h/3, k’= 3h in (17), and 
the stress is given by: 


= 3d 
y= Hoa [((v—h) log r,—3(a— $h) log r,, 
+3(a+4h) log r_,,—(@+h)logr_,]. . . . ~ (20) 


This has the maximum value 3-3ud/7(1—oc)h compression at x=0 and 
changes to tension at about a=0-47h. 
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are smaller as for k=0-2h, k'=0-1h the maximum 


If both parameters zz 
Sud/r(1—o)h and yy changes 


compressive stress increases slightly to 3- 


sign at about 2—0-6h. 
If & is zero as in (18) the stress 7y becomes infinite at +h, but remains 


finite at e—0 with the value 3-4ud/7(1—o)h and changes sign at about 
0-7h. 

If & is zero and k’ —>h in (18) the graph of v’ becomes a smooth curve 
with no straight segments. The dislocations are situated only at x=+h, 
with infinite tensile stress at these points. The compressive stress at 
a=0 ig 2ud/m(1—o)h, and 7 changes sign at x=0-75h. 

If the relative displacement of atom centres is required, then elastic 
terms must be added to v'(x). If the inter-atomic distance normal to the 
slip plane is a, then the relative displacement Av is calculated from the 


equation. 


Vo=v(x, 5a) u(z, 5a) 


and differs from v'(a) by a term depending on a. 


§ 5. Strarn EnerGy anp DENSITY CHANGE 


The strain energy per unit length in the z direction of a large elastic 
body containing the dilated slit is calculated from the surface integral 
in the usual way. But the field of each slit considered is that of the stress 
function Clog r+-cos? @) at large distances. Such stresses and displace- 
ments decrease so rapidly for increasing r that their contribution is 
negligible if the radius R of the outer surface is sufficiently large. The 
remaining term of the strain energy is given by: 


rh 
=—| Gj.v'dx SS coe. Cet 


where v’ is given by a chosen form of (16), and yy is deduced from (6) 
and (17). 
_ This quantity # may be calculated for any chosen values of k and k’. 
It may even be calculated for zero k since the product v’ . 7 remains 
finite as 4% approaches an infinite value at x=--h. For k zero, the 
strain energy is found to have the value 0-74yd2/7(1—<) for k’=0-1h, and 
the larger value pd?/7(1—o) if k’=h. For evenly spaced dislocations 
with k=$h, k’=th, the strain energy is 0-8ud?2/7(1—o), while for smaller 
parameters k=0-2h, k’=0-1h, the strain energy is. sligl 
peer eset gy is slightly smaller, 
When the field of the dilated slit is used to represent local expansion 
at an edge dislocation it is found that the strain energies of slit and 
dislocation may be simply added. This is because there are no Cross 
terms in the strain energy of the combined system, one having relative 
displacement w’ and zero yy on the x axis, the other a non Hatt v’ but 
vanishing x7. 
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The change in density of an elastic body due to the presence of an 
edge dislocation is zero according to first order elasticity theory. But 
if the lattice is dilated at the centre by atomic misfit, even first order 
elasticity admits a change of volume, and therefore of density. 

The increase in volume VV per unit length in the z direction of all 
the elastic material initially bounded by a cylinder of radius R is found 
by integrating the radial displacement round the outer surface : 


v= | chido whl Mae te ad, ~ Lwfo0N 


where w,=w cos 6+-v sin 6 
C 
2uR 
by eqn. (10). This gives 


[1—2(1—o) cos 26] . . . . . (23) 


VV =rC/p 
where the constant ( is determined by the form of the dilation as shown 
in §4: 
= p.d(3h?— 3hk+k?—k’?) 
Sn o)\2n—-’ 


But the u, of eqn. (10) applies to a body which extends to infinity, 
so that certain stresses act at the surface of radius R. If this surface 
is supposed free, it is necessary to add a stress function U’, as follows : 

ij C 2 qa 6 1 2 ie 
Urn= pe r? cos 20+-3r2— Re 08 20 


which contributes a further radial displacement 


C' f[1—2o 
v= 5 —(1 a) eos 20). 


On adding this to (23) the change in volume is found to be 


VV =27C(1—o)|p. Se gee, on. Gl 4) 


This VV. as might be expected, is equal to the volume of imaginary 
material added within the slit. 
If the average density of the body is p, the change in density is given 


approximately by 
VV 
ee Pe 


= — 2pC'(1—o) wR. 


§ 6. CONCLUSION 
It has been shown how to calculate the elastic stress field surrounding 
a slit dilated by an arbitrary distribution of enclosed material. At a 
distance the field is that of a particular type of centre of pressure, the 


P.M, B 
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actual values of the stresses depending on the contours of the inclusion. 
For a smoothly tapered inclusion the stresses remain finite everywhere, 
but an angular form causes singularities at isolated points. 

The strain energy of a large body containing such a slit may be 
calculated with accuracy, but the density change, as calculated in § 5, 
is less reliable. As pointed out by Stehle and Seeger (1956) the second 
order terms of the elastic dilation are no longer negligible when integrated 
over a large volume, and their total may even exceed those of first order 
theory. This effect makes the displacements of the present paper 
unreliable at large distances, but does not affect the components of 
stress or strain. 

It seems clear that the stress field of an edge dislocation in a crystal 
lattice must include terms of the type discussed here. At large distances 
such terms may often be negligible, but they may have considerable 
influence on the migration of solute atoms and the formation of precipitates 
on the slip planes. 
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ABSTRACT 
Earlier results obtained at Bristol from a study of nuclear interactions 
produced by heavy nuclei of the primary cosmic radiation passing through 
nuclear emulsions have been amplified and confirmed. Revised values are 
given for the fragmentation probabilities and mean free paths which are of 
greater statistical weight and in good agreement with those reported 
previously. Values are derived from the fragmentation probabilities in air 
and in hydrogen which probably represent upper limits to the true values. 

The implications of the results are discussed. 


$1. IntTRODUCTION 


THE experiment reported here is a continuation of a previous paper by 
Fowler, Hillier and Waddington (1957), which will be referred to as 
Paper I. In that paper the results were given of an analysis of 317 
nuclear interactions produced by heavy cosmic ray nuclei passing through 
nuclear emulsions exposed at high altitude. The present paper reports the 
result of a similar analysis made on a further 470 interactions. 

As before, the incoming nuclei have been separated into charge groups 
distinguished in the following manner : If Z is the charge of a nucleus, then 
L nuclei have 3<Z<5; WM nuclei 6<Z<9; A’ nuclei 10<Z<19; 
VH nuclei Z>20 and H nuclei Z>10; so that H nuclei include H’ nuclei 
and VH nuclei. Thus, in order to determine the fragmentation proba- 
bilities, i.e. the probability that a nucleus of a given charge group 
should produce, in a nuclear interaction, a secondary nucleus of the same 
or lighter charge groups. it is necessary to determine only the charge groups 
to which these nuclei belong. In practice it is generally possible to measure 
the charge of a nucleus to within one charge and thus the errors introduced 
by incorrect charge identification are usually small. For example. since 
the completion of Paper I further charge measurements have been made 
on many of the 1 and M nuclei observed in that experiment (Waddington 
1957 a), which show that less than 19%, of the carbon and boron nuclei 
were assigned to the incorrect charge groups and that the error in the 
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relative numbers of 1 and M nuclei was less than 4%. The charge 
determinations on the L and M nuclei in the present experiment should be 
more accurate than those in Paper I, since they were made by an observer 
with greater experience of making these rather subjective measurements. 


Fig. | 
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Flight curves of F and J stacks. Altitudes expressed in thousands of feet. 


From the combined data of the two experiments an attempt has been 
made to determine upper limits to the fragmentation probabilities in air 
and in hydrogen. The values obtained strongly support the view that 
the observation of a ratio of L to M nuclei greater than about 0-25 at 
10 to 15 g/cm? of residual atmosphere is incompatible with the assumption 
that there are no / nuclei incident at the top of the atmosphere. Further- 
more, it appears that these values of the fragmentation probabilities in air 
can only be raised by reducing those in hydrogen and thus increasing the 
significance of any observed primary 1 nuclei. 


§2. EXPERIMENTAL PROCEDURE 
Of the 470 interactions found in the present investigation, 84 were 
observed in the same stack of nuclear emulsions as those analysed in 
Paper I, the # stack, while the remaining 386 were observed in another 
stack, the / stack. exposed under closely similar conditions. Details of 
the exposures and physical dimensions of these stacks are given in table 1, 


while the flight curves are shown in fig. 1. Both these exposures were 
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made over Northern Italy, where the geomagnetic cut-off energy is about 
1-5 Bev per nucleon (Fowler and Waddington 1956), and as a result the 
great majority of the nuclei entering the emulsions were moving with 
relativistic velocities. 


Fig. 2 
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Charge spectrum of VM and L nuclei observed in J stack. Secondary fragments 
are shown shaded. 


Emulsions of the / stack were scanned, along a line parallel to the top 
edge, for the tracks of multiply charged nuclei. These tracks were then 
followed through the emulsions until the particle interacted or left the 
stack. Any fast secondary fragments heavier than «-particles emitted 
from these interactions were also followed until they in turn interacted 
or left the stack. In the # stack, where the interactions produced by 
nuclei entering the top edges had already been analysed in Paper I, a scan 
was made parallel to the side edges and of the end plates. These coum 
were made principally for H nuclei. Throughout, tracks were only 
accepted if they had a projected length of greater than 4mm per emulsi 1 
and a zenith angle of less than 60°. : ais 

The charges of the L and M nuclei were determined by 3-ray counting 
with similar procedures to those employed previously Caddie 1957 a 
The results of these measurements in the J stack are shown in fig. 2 aa 
charge assignments shown in this figure were based on the a eecun of 
charge-indicating interactions. Similar measurements were made on the 
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appropriate tracks in the H stack and the results are shown in fig. 3. 
In this case the charge assignments have been calculated from the 
following relation determined previously for this stack (Waddington 
1957a): Z=3-76 (N,—0-46)1, where N, was the 5-ray density per 
100 uw. 

The charges of the majority of the H’ and VH nuclei have been deter- 
mined with a photo-densitometer, and details of the procedures employed 
will be published in a forthcoming paper by Hillier and Rajopadhye. It 
is estimated that these determinations have a probable error of half a 
charge, if it is assumed that the highest charge peak with an appreciable 
abundance observed in their charge spectrum is due to iron nuclei. Ina 
number of cases it proved impracticable to make these densitometer 
measurements and the charges have been estimated from a 5-ray count. 
These latter determinations are not claimed to be very accurate, but do 
generally enable the separation of the nuclei into different charge groups. 
Particles identified in this manner are indicated in the table of interactions 
(table 2) by asterisks. 


NUMBER 
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Charge spectrum of J/ and / nuclei observed in # stack. Secondary fragments 
are shown shaded. 


All the charge determinations made in primary nuclei were made near 
the point at which they were detected, while all those on the secondary 
nuclei were made near their parent interactions. In a few cases, less than 
ten, this procedure resulted in an obvious anomaly between the measured 
charge and that apparent near the interaction. These anomalies, when 
confirmed, were found to be due to the presence of undetected interactions, 
which were generally either of the type in which only particles with mini- 
mum ionization were produced together with the heavy fragment of 
reduced charge, or to interactions which occurred very near to an interface 
of the emulsions. The occurrence of these events indicated that not all 
interactions of these types can have been observed. However, it appears 
very unlikely that an appreciable proportion were missed, because of the 
small number found in this manner. This conclusion is supported by the 
close agreement of the observed mean free paths with those found by 
other workers (Noon and Kaplon 1955). While no attempt has been made 
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to improve the values quoted in Paper I for the mean free paths of L and M 
nuclei, new values of 39-6+5:4 g/em? and 26-8-+-5:5 g/cm? have been 
obtained for the H’ and VH nuclei. These values, taken in conjunction 
with those quoted in Paper I}, give 38-7--3:7 g/em? and 31-24 4-4 g/cm? 
respectively for the mean free paths. 

Every interaction observed was classified in terms of the number of 
black and grey prong emitted, V,, the type of secondary fragment, if any, 
the number of secondary «-particles, n,, and the number of fast singly 
charged particles, »,. In Paper [it was shown that for interactions of the 
type observed in this experiment none of the secondary fragments, LS2, 
were ejected with angles of greater than 10° to the line of flight of the 
incoming nucleus. For this reason all the tracks of heavily ionizing par- 
ticles emitted within this angle were closely studied. and if necessary traced 
through the stack, to determine whether they were due to secondary 
fragments or to slow evaporation particles from the target nucleus. It 
is, therefore, very unlikely that any secondary fragments can have been 
missed in this work. 


§ 3. RESULTS AND ANALYSIS 

The detailed features of all the interactions observed in this experiment 
are shown in table 2, with the exception of those of the type where one 
black prong is the only product. These interactions have been neglected 
in this work, as they were in Paper I, because of their low efficiency of 
detection. Anexamination of the data presented in this table showed that 
these results were consistent with those of Paper I. For this reason, in 
what follows, all the data have been combined. 

While a knowledge of the fragmentation probabilities in emulsion, P.,,,,. 
enables the observed flux values to be extrapolated to the top of the de- 
tecting emulsions, these are not the correct values to use for an extrapola- 
tion to the top of the atmosphere. At present there appear to be no 
experimental results on the appropriate fragmentation probabilities in air 
and very few on those in hydrogen. Previously the values in air have 
either been obtained from those in emulsions by calculation from a 
geometrical model of the nucleus (Noon and Kaplon 1955), or from the 
fragmentations produced in glass (Bradt and Peters 1950). Im this paper 
an attempt has been made to separate the interactions with the light 
elements in emulsion—carbon, nitrogen and oxygen—from those with the 
free hydrogen and from those with the heavy elements silver and bromine. 
This separation has been made on a statistical basis and depends on the 
number of slow particles emitted from the target nucleus, NV a 

The observed interactions were divided into three classes : 


(i) X-stars. These stars will be characterized by having V,—0 or 1. 
If there is a black or grey prong then for the star to be included in this class 
the prong must not be at an angle of greater than 90° to the line of flight 


+ The mean free path for VH nuclei was erroneously quote 


d in this paper, and 
should have been 35-1 6-8 g/em?, a : 
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of the incoming nucleus. If V,—0 then the sum of the charges of the 
particles produced in the star must be at least one charge greater than that 
of the primary nucleus and at least one of these particles must be singly 
charged. Stars of this class may be due either to interactions with free 
protons in emulsion, p-interactions, or to collisions with complex nuclei 
in which at most one slow charged particle is emitted from the target 
nucleus. 


(ii) R-stars. These are characterized by having N,<7 but do not 
include X-stars. They may be due either to interactions with light 
nuclei in emulsion, /-interactions, or to interactions with heavy nuclei 
which do not seriously disrupt the target nucleus. 


(iii) 7’-stars. These are characterized by having V,>7 and are due 
to interactions with the heavy nuclei in emulsion, h-interactions. 

Before it is possible to attempt a separation of the stars in these empirical 
groups into those due to p-, /- and h-interactions, itis necessary to determine 
the true proportions of such stars which will be produced by a given 
incoming nucleus. To do this for collisions betwene complex nuclei the 
interaction cross sections, o, have to be calculated for the different 
constituents of the emulsion. A model which permits this calculation has 
been proposed by Bradt and Peters (1950) to fit the observed mean free 
paths of different nuclei in glass and brass. These authors proposed an 
‘overlap ’ model such that o of anucleus with a radius R , when approached 
by a nucleus with radius R, was given by : 

Wie fiw al ee ee ee oe 8) 
and AR had the empirical value of 0:85x10°-8 cm. Now R=r,A1? 
where A is the atomic weight, and in 1950 7, was assumed to be about 
1:45x10-8 cm. This relation has now been checked and found to be 
applicable over a wide range of media and incident nuclei, ranging from 
x-particles in air (McDonald 1956) to H nuclei in lead (Hisenberg 1954). 
If the presently accepted value of 7°) (1-20 x 10-18 cm) is used instead, then 
AR must be varied to fit all the data, decreasing with increasingly heavy 
incident nuclei, as might be expected from such a model. Thus the general 
applicability of eqn. (1) must be regarded as fortuitous. 

While eqn. (1) appears to be valid for all incident nuclei between helium 
and iron and all target media between air and lead, it is not valid for 
incident protons, nor presumably for a hydrogen medium. However. 
values of o for 860 Mev protons in various elements have been measured 
by Chen ef al. (1955) using an attenuation method, while o for proton— 
proton interactions at the same energy is known from the work of 
Shapiro et al. (1954). These results are shown in fig. 4, which shows that 
the experimental values lie on a smooth curve. This curve has been used 
to obtain the cross sections of the various constituents of emulsion, and a 
mean free path for protons in emulsion of 37 cm calculated. This valueisin 
close agreement with that found experimentally and may be contrasted 
with the value of 28 cm found from eqn. (1). 
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hus. when calculating the partial cross sections in emulsion for p-, 
l- and h-interactions, eqn. (1) was used for the complex target nuclei, - 
while the data of fig. 4 were used for the proton component. Due to the 
small proportion of hydrogen in emulsion this procedure hardly affects the 
total cross section in emulsion. These proportions for different incident 
nuclei, and the total mean free paths, are shown in table 3, together with 
the observed proportions of X- R- and 7'-stars. 

From the data shown in these tables it is possible to determine the 
number of /- and h-interactions in the X- and R-stars—all the 7’-stars 
being due to /-interactions—provided that it is assumed that those X-stars 
which are not p-interactions have the same proportion of /- and h-inter- 
actions as the entire sample. Fragmentation probabilities, P',,, were 
then calculated for the 7'-stars and the further assumption made that these 
values were those for h-interactions, P” 
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The inelastic cross section, o, as a function of the cube root of the atomic weight 


for protons on various elements (from data of Chen et al. (1955 
Shapiro et al. (1954)). pt a. (1955) and 


Since the 7’-stars are all interactions in which the heavy target nucleus 
has been considerably broken up, and must, therefore, be generally due 
to moderately central collisions having small impact parameters the 
incoming nucleus must also generally have been severely broken up. Ke a 
result, the probability of the incoming nucleus leaving a residual nucleus, or 
secondary fragment, should be less for these 7'-stars than for those Dati. 
actions, with their larger impact parameter, included among the R- and 


X-stars. Thus, the values for P”,,, calculated from the T-stars would 
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appear to be underestimates of the true probability that a secondary 
fragment will be produced in a /-interaction., 

Fragmentation probabilities. assumed to be those for /-interactions. 
P',,, were then obtained from the R-stars by making a correction for that 
proportion of them which were due to h-interactions by using the values of 
P', found above. A similar procedure for the X-stars, correcting both 
for the /- and h-interactions, then gave the fragmentation probabilities for 
p-interactions, P?,,. These values of P!, and P”, will be upper limits 


mine weet 


to the true values. This is true for P’ because the values of P”. used 


mere in 
to correct the data were lower limits to the true values, and for P”,,, because 
the main correction is also due to P’,, although in this case there is 
admittedly also a correction due to P',,,, which is itself an upper limit. 

This derivation is dependent upon the assumption that the fragmenta- 
tion probabilities derived from the 7'-stars are always less or equal to those 
of the other h-interactions among the R- and X-stars. It is not possible 
to prove that this assumption is correct, although it is very plausible. 
However. for the most important fragmentation probabilities, those into 
M and L nuclei, the correction introduced by using P”,, is small. For 
example, P’, has its greatest value for P!,,—0-15, yet even if it were 
assumed that P3,,=0 for R-stars, P),, would only be raised from 0-21 to 
0-25, an increase much less than the statistical uncertainty. The implicit 
assumption has also been made that the fragmentation probabilities for 
/-interactions among the X-stars will not be very much greater than for 
those among the #-stars. Since the proportion of /-interactions among the 
X-stars is, in every case, less than 10°, of those among the #-stars, unless 
this assumption is very seriously in error the values of P!,,, should not be 
affected. Moreover, it should be noted that any variaion in P’,,, must 
affect P?,, in the reverse manner. For example, the results of bombarding 
carbon with 90 Mev neutrons (Kellogg 1953) suggests that P/),,~1, 
implying that the quoted values of P’,, might be too small. If this is 
correct then the quoted values of P4,,, must be too high. 

The values obtained for P!,, and P”,,, are upper limits to the fragmenta- 
tion probabilities in air and in hydrogen, and can be used to extrapolate 
observed flux values through these media. The majority of expriments on 
the relative abundances of 1 and M nuclei have used detectors exposed 
between 5 and 20 g/cm? of air below the top of the atmosphere. Hence the 
effect of passage through air on the ratio of 1 and VM nuclei has been calcu- 
lated for two cases: (a) when the flux of / nuclei at the top of the 
atmosphere, J,., is zero, and (b) when -/;. is 0-4 of the flux of M nuclei, 
J,,-. This calculation was made using the simplified diffusion equations 
reported previously (Waddington 1957 b), which neglect the effect of 
consecutive interactions, and using the interaction mean free paths shown 
in table 3. The results of these calculations are shown in fig. 5, together 
with the experimental point found by (Waddington 1957 a) under 12 g/cm? 
of residual atmosphere. It can be seen that while the statistical errors 
are such that at depths greater than 20 g/cm? itis difficult to distinguish 
between cases (a) and ()), for lesser depths there is a marked divergence, 
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It is clear that the experimental point shown is not consistent with 
J 7e/J y=, particularly when it is remembered that this extrapolation 
probably represents an overestimate of the correction for secondary 
production. 


Fig. 5 


TO M-NUCLE! 


L-NUCLE| 


RATIO OF 


The ratio of to M nuclei as a function of the depth of residual atmosphere in 
g/em?, with the limits due to the statistical errors of the derived values of 
the fragmentation probabilities. 


When extrapolating the fluxes at the top of the atmosphere through 
interstellar hydrogen it is necessary to consider the effects of consecutive 
interactions and to use the diffusion equations published by Kaplon et al. 
(1954). Unfortunately the statistical weights of the P?,,, values are so low 
that only tentative conclusions can be drawn regarding the propagation 
of these nuclei. However, it does appear that, assuming J7./J ,.=0-4, 
the observed charge spectrum is consistent either with a pure H nuclei 
source after passing through about 20 g/cm? of hydrogen, or a source with 
no L nuclei after passing through about 8 g/cm? of hydrogen. If a density 
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of 10-24 g/em3 is assumed for the interstellar hvdrogen, i.e. the cosmic rays 
are at least confined to the galactic disc. then these two cases correspond 
to lifetimes of 20 million and 8 million vears respectively, which are con- 
siderably longer than the values currently accepted (Rossi 1955). If 
cosmic ray particles are confined within the galactic halo then these life- 
times will be as much as two orders of magnitude greater. It should be 
noticed that if J,./-J,). is in fact less than 0-40 because the values eid gee 
obtained in this paper are too low, then the quoted values of P%,,, must be 
too high and a smaller J ,/J ,,. ratio will still imply an appreciable lifetime. 

Detailed values of the fragmentation probabilities in hydrogen are of 
considerable importance in our understanding of the significance of the 
observed cosmic ray charge spectrum. It is to be hoped that the 
theoretically quite simple, though experimentally difficult, experiments 
to find these parameters will soon be undertaken. 


§ 4. CONCLUSION 


The most important conclusion that we reach from this work is that the 
several recent experiments, reported at the Varenna Conference, which 
give J,/J,, 0-4 at 8-15 g/cm? of residual atmosphere are not compatible 
with absence of L nuclei at the top of the atmosphere. 
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ABSTRACT 


Negative K-mesons captured by the heavy nuclei in nuclear emulsion form 
K-mesie atoms. It is shown that absorption occurs in the nuclear surface 
and takes place when the K~-meson is in a 5g state. Observations on capture 
events producing a visible 2~-hyperon, and on the ratio of 7+ to 7~-mesons 
created can give information on the relative neutron and proton densities 
in the part of the nuclear surface where the total density is less than 10% 
of the density at the centre of the nucleus. Experimental evidence suggests 
that in this region of the nucleus the density of neutrons is approximately 
equal to the density of protons. 


§ 1. IyTRODUCTION 


JOHNSON AND TELLER (1954) have suggested that as a result of the 
modification of the proton charge independent potential by the Coulomb 
potential, the density of neutrons is considerably greater than the density 
of protons at the surface of a heavy nucleus. The radial density distri- 
butions of neutrons and protons in heavy nuclei have been compared 
experimentally by Hess and Moyer (1956) who studied the deuterons 
produced in (p,d) and (n, d) reactions at large angles to the incident 
beam. Using 300 mev nucleons they found evidence for a possible 
excess of neutrons in the surface regions of heavy nuclei. For lead this 
excess could take the form of a surface layer of neutrons 0-8 x 10-3 cm 
thick. On the other hand, Abashian et al. (1956), by measurements 
on the absorption cross sections of 700 Mev 7+ and 7 -mesons in lead 
find no evidence for a neutron excess in the nuclear surface. These 
methods however, are not sensitive to the relative numbers of neutrons 
and protons in the extreme part of the nuclear surface where the total 
density is less than 10% of the density at the centre of the nucleus. 
The purpose of this note is to point out that this part of the nuclear 
surface can be investigated by observations on the absorption in heavy 
nuclei of slow negative K-mesons. 

The K--meson is captured in a Bohr orbit of very high principal 
quantum number and then makes electromagnetic transitions to states 
of lower energy. The theory of Wightman (Bethe and de Hoffmann 1955) 


+ Communicated by the Author. 
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shows that absorption by the nucleus takes place from a state with the 
highest possible angular momentum. In § 3 it is shown for a nucleus 
with A—94 and Z—41, that this is the 5g state. In the Schrodinger 
equation for the bound K~-meson the nucleus is represented by a complex 
potential V=U+iW. Unless U is attractive and |U|> |W |, the 
K--meson wave function for a 5g state must have extremely small 
values inside the nucleus and the absorption must take place on the 
surface. In §2 the range of values that U can have is investigated 
semi-quantitatively. 

In almost all cases the K~-meson is absorbed by a single nucleon 
producing a hyperon and a z-meson (Dilworth 1957). If the inter- 
action is charge independent, the amplitudes a, and a, for the production 
of a 2-hyperon in states with isotopic spin 7=0 or | are related by 
a,—=a,(0-8e7) for the absorption of slow K~-mesons by free protons 
(Alvarez 1957). The 7’'=1 amplitude, b,, for production of a A°-hyperon 
satisfies |b, |?=0-07 | a, [?. The normalized transition probabilities for 
possible reactions, calculated using these relations, are given in table 1. 


Table 1. The normalized transition probabilities for the different modes 
of absorption of a slow K~-meson by a single free proton or neutron, 
calculated from the experimental relations found by Alvarez 
(1957), are shown 


; Normalized transition 
Reaction salen 
K-+p > tt at 0-14 
QT 0-28 
2M -|-479 0-1 i! 
A%17,0 0-02 
Si Ue 0-20 
2 +79 0-20 


Ae a= 0-05 


The transition probabilities for absorption in nuclear matter may be 
different because the relations between the amplitudes Go, Gy and. 5 
depend on the relative momentum and angular momentum of ‘ine 
K™-meson and the nucleon. This dependence is unknown but is assumed 
to be weak. In support of this it should be mentioned that the relative 
transition probabilities observed for the absorption of fast K--mesons 
by free protons (Rosenfeld 1957) agree with the results of table 1 within 
the present large experimental errors. Final state interactions other 
than Coulomb, should not affect the transition probabilities if the Bbc 


oe takes place on the surface of the nucleus where the nucleon density 
is low. ; [ 
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Using the results of table 1, the expected ratio of the number of 
7* to the number of = -mesons produced for absorption in a heavy 
nucleus with 4=94 and Z=41 and neutron and proton radial density 
distributions equal, is R=0-61. In addition, the ratio of the number of 
+ -hyperons produced by the reaction 


K-+n-> 2--+7° or a pce ere en rae Bie 8 
to the number produced by the reaction 
Ree Sell eeu ne ey eats (Lie) 


should be S=0-94. Experimental values exist for both ratios (Dil- 
worth 1957). In §3 they are shown to depend strongly on the relative 
numbers of neutrons and protons in the nuclear surface. In § 4 the 
values of R and S calculated as a function of the neutron excess in the 
nuclear surface are compared with the experimental values. 


§2. THe K--Nucievus PoTENnTIAL 


In this section an estimate is made of U, the real part of the K~-nucleus. 
potential. As an approximation to this we consider the energy level 
shift of a K--meson interacting with an infinite nucleus. Coulomb 
interactions are neglected. The absorption of a K~-meson by two 
correlated nucleons is experimentally observed to have low probability. 
One reason for this is that in the centre of mass system of the two 
nucleons, each nucleon would need to have a momentum of 460 mev/c 
to produce a real final state of a nucleon and a hyperon. The prob- 
ability of such high momentum is small. However, elastic scattering 
through an intermediate state of a nucleon and a hyperon can still make 
an important contribution to U. It will be shown that this contribution 
is always positive in sign. The approximations to be used here are 
almost identical with those used by Deser et al. (1954) and by Brueckner 
(1955) to estimate the energy level shift of the s-state in 7-mesic atoms. 

The energy level shift of a single meson interacting with a system of 
nucleons has been calculated by Riesenfeld and Watson (1956) and 
Watson (1957). A first approximation to their result is 


U=U,+0, 

dy (BIRIg al RIB) 
(277)? wy +2H, —E,(¢)—E£(4) 
U,=N.Re<f | t(0) | 8). h=c—1, 


Here P denotes the principal value of the integral. & is an operator 
which destroys a K~-meson and two nucleons and creates a nucleon and 
a hyperon. The |¢) are a complete set of states of one nucleon and one 
hyperon, and |f) is the nuclear ground state. In the calculation of 
U,, the momentum transform of the two-nucleon wave function 4(q) 
used is the Gaussian form proposed by Brueckner (1955). In the present 


C2 


U;=N,P | 


36 P. B. Jones on the Absorption of Negative K-Mesons 


case this term gives a negative contribution to Us only for values of 
4(q) for which q is greater than 460 mev/c, and so the sign of U, is not 
affected by the momentum dependence of F and #(q). We reasonably 
assume R to bea constant. w, H and H, are the energies of the K~-meson, 
nucleon and hyperon respectively. N, is the effective density of corre- 
lated pairs of nucleons. The mean free path in nuclear matter for 
absorption of a K~-meson by two nucleons is A=(N 50) where o 1s 
the cross section. Two-nucleon processes involving the production of 
virtual K+-K~ pairs should be unimportant and will not be included. 
In this cage U,=45000(A)~! mev, for small values of v, the K~-meson 
velocity in the centre of mass system of the two nucleons. A is in units 
of 10-3 cm. If absorption by a pair of nucleons occurs in 5% of all 
cases, A=10-12. cm and setting v=0-1, U,;—45 Mev. 

To estimate U,, we relate the K~-nucleon scattering matrix at zero 
energy, 1(0), to the zero-energy s-state cross section for K~-p elastic 
scattering using for this an extrapolated value of 40 mb (Alles e¢ al. 1957, 
Rosenfeld 1957). Assuming that the real parts of the K-—p and K-—n 
zero-energy s-state scattering lengths are both positive and equal in 
magnitude, we find U,=-—70 mev for a K~-meson in nuclear matter of 
density N,—0-17 x 108° nucleons per cm’. This is the extreme negative 
value that can be contributed to U, unless the K-—n cross section happens 
to be greater than 40 mb. ; 

The actual value of U is not important for present purposes. However, 
these results combined suggest that U cannot be large and attractive, 
and therefore that absorption occurs in the nuclear surface only. In 
addition there is some experimental evidence that the real part of the 
optical model potential for K~-mesons is attractive and has about the 
value estimated here for U (Alles et al. 1957). 


§3. CALCULATION OF R AND S 


The ratio & of 7+ to 7~-mesons produced in the surface of the nucleus, 
and S, the ratio of 2~-hyperons produced by reaction (1.1) to 2~-hyperons 
produced by reaction (1.2) are calculated in this section as functions 
of the neutron excess in the nuclear surface. All calculations are for a 
heavy nucleus with A=94 and Z=41, which represents the elements 
Ag, Br and I in nuclear emulsion. The light nuclei are represented by a 
nucleus with Ad=13-8 and Z=6-9. 

In first order perturbation theory the mean life 7, for nuclear absorption 
of the A~-meson from a state with unperturbed wave function @ is 
given by (rt) *=<(¢|W|¢). W is the imaginary part of the complex 
K~-nucleus potential and is taken to be +100 Mev. This is obtained 
from the mean life in nuclear matter which is 0-6 x 10-23 seconds estimated 
using the measured absorption cross section, in hydrogen, for K~-mesons 
of momentum 10-170 mev/e (Rosenfeld 1957). It has been assumed 
here that the free particle absorption cross section is unaltered in nuclear 
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matter of normal density}. In table 2 the approximate mean lives of 
several states are given both for transitions (4/—-+1) to a state of lower 
energy with the emission of a y-ray (Dirac 1947) and for nuclear absorp- 
tion. It is seen that absorption takes place from the 5g state. It is 
worth mentioning that some absorption from the 4f or 6h states would 
in no way alter the results to be obtained here. 


Table 2. The approximate mean lives of bound states of captured 
K -mesons are given both for 4/=-+1 y-ray transitions to states: 
of lower energy and for nuclear absorption. The mean lives for 
nuclear absorption are calculated using a K~-meson mean life 
in nuclear matter of 0-6 x 10-5 seconds. The charge distribution 
found by Hofstadter (1956) is assumed equivalent to the nucleon 
density distribution 


x Mean life for electro- Mean life for 
State F os : 
magnetic transitions nuclear absorption 
6h 1-0 >< 1077 sec 10375 sec 
5g 39x 10-17 sec 0-4 x 10-1!’ sec 
4f 1-3 x 10-** sec 0-25 x 10-1* sec 


u(r)=r¥(r) has been found numerically for U-+->~W=0+1001. Yr) is 
the K--meson radial wave function for the 5g state. The nuclear 
potential has been taken as 

V=(U+iW)lexp {(r—a)/z,} +1 

where a=5:5<10-% cm. z,=0-55x10-%cm, the value found by 
Hofstadter (1956) from measurements on the nuclear charge distribution. 
The form of u(r) is not sensitive to the exact values of the parameters 
a, U and W. The behaviour of w(7) depends mainly on the high angular 


momentum. 
The probability of the absorption taking place at a distance r from 


the centre of the nucleus is 
P(r)= const. | u(r) |?{pnPn(r)+(@p/on) ppl p(7)}- 
pn and pp are the neutron and proton densities at r=0. wn and wp are 


the transition probabilities for absorption of the K~-meson by neutrons 
or protons (see table 1). The proton form factor (Hofstadter 1956) is 


Fy(r)=1 for 0<r<(c—Zs) 
(c+23—1)/2z5 (C—23) <r< (€+2s) 
0 (ctz,)<r 


+ For the optical model of the nucleus this approximation is good (Lane 
1957). Watson (1957) points out that the potential V is essentially the optical 


model potential. 
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Two neutron form factors have been used. The first has the same form 
as Fy(r) except that the parameter c is replaced by c+4c. Ac is then a 
measure of the neutron excess. For the second neutron form factor, 
Fy(r)=F p(r) for r<(c+z 3). However for (c+23+ Ah) > r > (c+2s); 
Fx(r)=0-05F p(0). This represents an outer shell, 4h in depth, of 
neutrons with extremely low density which would not be detectable by 
other methods. 

For the light nuclei the suggestion of Johnson and Teller, if correct, 
would lead to de~10-“ cm only, for a neutron potential of —35 Mev 
and the charge distribution found by Hofstadter. It is therefore assumed 
that for the light nuclei #y»=/'p. The ratio of the number of K”-mesons 
captured by light and heavy nuclei is taken to be 3/7. 

In the case of the ratio R, it is necessary to consider the interaction 
of the pions with the nucleus. The hypothesis of charge symmetry 
implies that 7+ and 7~-mesons are scattered and absorbed equally in 
self-conjugate nuclei. For the heavy nuclei this is not true. To estimate 
the correction we use the pion absorption model of Brueckner e¢ al. (1951), 
and the phenomenological theory of pion production by two nucleons 
of Rosenfeld (1954). The correction arises from absorption of the pion 
by two protons or two neutrons which are assumed to be in a singlet 
s-state. The difference of the mean absorption cross sections per nucleon 
for «+ and 7 -mesons is o(+)—o(—)=0-15I7?(1+-7?)! mb. 7 is the 
momentum of the pion in the centre of mass system of the two nucleons 
in units of wc (uw is the pion mass). For the correction we have used 
I'=4}, which results in a change of less than 4° in R. Geometrical 
corrections to k and S, present only if there is a neutron excess at the 
nuclear surface, are at the most a few per cent. 

The transition probabilities of table 1 can now be used to calculate 
Rand S. In figs. 1 to 4, R and S are shown as functions of dc and Jh. . 
In all figures the curves (a) show the values of R and S for capture by 
the heavy nuclei alone. The curves (b) show the values of R and S 
resulting when K--mesons are captured by both light and heavy nuclei 
in the ratio 3/7. 

§4, Discussion 


The calculated values of R and S are seen to depend strongly on the 
relative numbers of neutrons and protons in the nuclear surface. The 
dependence is of course greater for the case of absorption of K~-mesons 
by the heavy nuclei alone. 

The value of S (Dilworth 1957) found from observations on K--meson 
capture events in nuclear emulsion which produce a visible 2, -hyperon, 
is 0-6. With neutron and proton density distributions equal, the predicted 


+ Fis defined (Francis and Watson 1953) as the ratio P/P,, where P is the 
probability, in nuclear matter, of a neutron and a proton being separated by a 
distance d, and P, is the corresponding probability for the ground state of the 


deuteron. d=h/p, where p is the relative momentum that the two nucleons 
need to be able to absorb a pion. 
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value is S=0-87. Figure 3 indicates that 4c, for the absorption model 


used here, is probably less than 0-2 10-18 em. For the second neutron 
form factor the corresponding limit is 4h <0-5 x 10-" em. 


Fig. 1 


to ma mesons 


+ 


Ratio of 


0:0 0:2 0-4 0-6 
Ac in units of 10°? cm 


The ratio R, of 7* to 7~-mesons produced in the nuclear surface is shown as a 
function of Ac, the neutron excess for the first neutron form factor. 
The curve (a) is for capture by the heavy nuclei alone, and curve (b) 
is for capture of the K~-meson by both light and heavy nuclei in the 


ratio 3/7. 
Fig. 2 
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R is shown as a function of 4h, the neutron excess for the second neutron 
form factor. 


R=0-64 is predicted when 4c=Ah=0. The true experimental value 
of Ris unknown. The existing value of R=0-20 has been obtained from 
the low energy end of the pion spectrum only, and so is very heavily 
biassed in favour of the 7~-meson (Dilworth 1957). 
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Fig. 3 


ORO) (ORe O4 0-6 


Ac in units of 10 cm 


The ratio S, of 2--hyperons produced by reaction (1.1) to 2~-hyperons produced 


by reaction (1.2) is shown as a function of 4c. 


Fig. 4 


HO) {(O¥4"  Oys (2 ie) 
Ah in units of 10 2m 


S is shown as a function of Jh. 
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An optical model may be used for the pion—nucleus interaction. The 
value of the ratio J” is very poorly known and so we do not take into 
account refraction and pion—nucleon scattering. The importance of 
scattering for pions of these energies is in any case much reduced by the 
exclusion principle. Using ['=4, we find for the K~-meson absorption 
model considered here, that the probability of escape of a charged pion 
is 0-66. Considering the large uncertainties in I" this compares well 
with the value of 0-53 found by Dilworth (1957). 

Finally, the results of table 1 deserve some comment. Alvarez (1957) 
has pointed out that the relation between a, and a, is based on the assump- 
tion that the absorption of a slow K~-meson by a free proton all takes 
place from one atomic state. Observations on absorption in deuterium 
should decide whether or not this assumption is good. In this case, 
absorption of the K~-meson by both nucleons should not, for the reasons 
mentioned in § 2, be important. 
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ABSTRACT 


The heat capacity of 160 g of diamonds has been measured in the temper- 
ature range 12-8° to 277°x with estimated accuracies of +20% at 13°K, 
+6% at 20°x, +0-8% at 100° and +0-2% for T>200°x. Fitting of 
the results in the region 7’< @p/20 to the low temperature expansion for 
the heat capacity gives @,—2219+20°K; this agrees excellently with 
@ (elastic) —=2240+-5°x, derived from elastic constants determined by 
McSkimin and Bond. Down to 75°x the results are in reasonable agreement 
with earlier measurements of DeSorbo, but at lower temperatures an 
anomaly reported by DeSorbo has not been confirmed. 


$1. [INTRODUCTION 


ALTHOUGH thermal, elastic and optical properties of diamond have been 
studied extensively for a number of years, rather large discrepancies in 
existing data are apparent. For example, the heat capacity of diamond 
has been measured below 300°K by Pitzer (1938) and by DeSorbo (1953), 
but the two sets of results differ systematically by up to 15°%. Large 
differences also appear in elastic constant data. A lattice frequency 
spectrum has been derived by Smith (1948), but its reliability can hardly 
be judged with the present experimental information. It was therefore 
thought desirable to make the new heat capacity measurements which 
are described in this paper. Since the first draft of the paper was 
completed, Burk and Friedberg (1957) have made a preliminary report 
on similar measurements in the region 7'<200°K. 

At temperatures less than about @,/20 the heat capacity of a solid 
is given by 


CSOT OT pel Oe ee ea 


The coefficients a, b and ¢ may be determined most conveniently from 
a plot of C/T? against T?. The coefficient a@ may also be deduced from 
the elastic constants of the crystal, and this provides a useful correlation 
between thermal and elastic properties. Theory and a number of 
experiments suggest that the coefficients b and c are usually positive 
for simple crystals (Blackman 1937, Barron and Morrison 1957), but 
DeSorbo’s (1953) low temperature results show @, increasing with 
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temperature just above 7'=0°x, thus indicating a negative value of b for 
diamond. This result, if correct, would be of considerable interest, 
especially as existing heat capacity data for other diamond structures 
(silicon, germanium and grey tin) suggest that b is small and positive 
for these elements. 

The present results show no maximum in @, at low temperatures 
and indicate that both 6 and ¢ are positive. The derived value of ©, 
agrees well with the value of O(elastic) computed from elastic constants 
recently reported by McSkimin and Bond (1957). 


§2. EXPERIMENTAL 


The calorimetric sample was approximately 160 ¢ of clear industrial 
diamonds ; their average dimension was about 3mm. Most of the 
diamonds were lightly coloured either yellow, brown or green, but none 
contained occlusions. Whether they were of type I or type II is not 
known ; this was not thought to be important because Robertson et al. 
(1936) had not detected a significant difference between the heat capacities 
of the two types. 

The measurements were made with a calorimeter assembly which had 
been used previously for experiments with alkali halide crystals. Details 
concerning the temperature scales, the handling of samples, etc., are 
given elsewhere (Berg and Morrison 1957). For the diamond experiments 
the calorimeter vessel contained 3-5 10~® moles of helium to promote 
heat exchange. This amount of helium was such that its desorption 
during a single measurement would not significantly affect the heat 
capacity even at the lowest temperatures. 


§3. RESULTS 


The original experimental results are recorded in the Appendix. The 
units of C, are cal/g atm deg, where | calorie has been taken as 4-184 abs 
joules. The apparent Debye characteristic temperatures have been 
calculated directly from C,,, since the difference C,—C,, is entirely 
negligible throughout the temperature region considered. A graph of 
apparent @,, against temperature is shown in fig. 1; the results obtained 
by Pitzer (1938) and DeSorbo (1953) are included for comparison. 

An estimate of the accuracy of the measurements can be made from 
previous extensive measurements on alkali halide crystals with the same 
calorimeter assembly (Berg and Morrison 1957). There is every reason 
to expect that the accuracy of measuring the total heat capacity of the 
assembly has remained the same, so that it is only necessary to take 
account of the much smaller contribution of the diamonds to the total 
heat capacity at temperatures below 200°K. In this way we obtain 
estimated accuracies for the measured heat capacities of +20° at 13°, 
+6% at 20°, +0-:8% at 100° and +0-2% for T>200°K. The results 
are consistent with a smooth 0, —T relation well within these estimated 
limits except for two points at the very lowest temperatures which are 
just outside them (fig. 1). 
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Fig. | 


7 ' 
0 50 100 150 200 250 300 
T (°K) 
The apparent Debye characteristic temperature of diamond as a function 
of the absolute temperature. 


O—present results. 
(1—Pitzer (1938). 
A—DeSorbo (1953). 
®—Smith (1948). 


$4. Discussion OF RESULTS 
4.1. Comparison with Earlier Measurements 


The present results agree rather well with DeSorbo’s down to 75°k. Their 
greater internal consistency is in part due to the larger specimen used— 
twice that used by DeSorbo (1953). Both sets of results show a con- 
sistently lower heat capacity than found by Pitzer (1938). The 
anomaly which appears in DeSorbo’s results at lower temperatures has 
not been confirmed; below 75°K, @) increases regularly and _ finally 
appears to level-off. 

DeSorbo’s anomaly may be due to his use of less perfect diamonds, 
but it could be explained more easily by a small systematic errer in the- 
heat capacity of the calorimeter vessel. The results would be very 
sensitive to such an error; for example, in the present experiments an 
uncertainty of 0-1°, in the heat capacity of the calorimeter vessel corre- 
sponded to an uncertainty of 3-3°% in the heat capacity of the diamond 
at 50°K. This is a more likely source of error than an undetermined 


effect of desorption of exchange gas, especially since DeSorbo obtained’ 
similar results with and without exchange gas. 
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4.2. Thermodynamic Functions 


The thermodynamic functions of diamond at higher temperatures 
are not greatly affected by the differences in heat capacity between the 
present results and DeSorbo’s. Values of S° and H°-H,° at 100°K for 
both sets of results are given in table 1. Higher temperature values 
may be obtained from DeSorbo’s table II after correction for the differences 
shown here in table 1. 


Table 1. Values of the Functions S° and H°—H,° for Diamond at 
100°K 


S° H°—H, 
(cal/g atm deg) (cal/g atm) 


De Sorbo (1953) 0-0192 1-3095 
Present experiments 0-0172 1-3211 


Cv/t3x10® cal./gm. atom deg.’ 


T2 x 1o73 


A graph of C,/T? vs T? for diamond. 


4.3. The Low Temperature Expansion for C, 


@,, the Debye characteristic temperature at 0°K, can be deduced most 
accurately from a plot of C,/T* against 7? (Barron and Morrison 1957). 
Such a plot for the diamond results is shown in fig. 2. The solid curve 
represents the best fit to the low temperature expansion for the heat 


capacity (eqn. (1)). The derived values of a, 6 and ¢ are given in the 
middle line of table 2. The other two lines of the table contain maximum 
and minimum estimates of the coefficients. The values of b and c depend 
critically on the value chosen for @, so that their magnitudes cannot 
be obtained very accurately. 
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Table 2. Coefficients of the Low Temperature Expansion for the heat 


Capacity 
a x 108 bx 101s Cx 1()*7 
(cal/g atm deg*) (cal/g atm deg®) (cal/gm atm deg’) 
4:35 5:5 11 
4-25 8:8 8 
15 13-0 4 


The form of fig. 2 proves to be much like that found for alkali halide 
erystals (Barron and Morrison 1957). Detectable changes in @, occur 
down to temperatures well below @,/50. The dependence of C, upon 
T up to T=@,/50 is given rather closely by the first two terms of eqn. (1). 


4.4. Comparison with Elastic Data 

The value of @, thus derived from the low temperature heat capacities 
(2219--20°K) provides a partial criterion for judging the accuracy of 
existing elastic constant data: these should yield a characteristic 
temperature, @(elastic), in agreement with @,. The most recent deter- 
minations of the elastic constants of diamond have been made by 
McSkimin and Bond (1957), and their results differ markedly from those 
of Bhagavantam and Bhimasenacher (1946) and of Prince and Wooster 
(1953), particularly for c,, and c,4,. Q@(elastic) is sufficiently sensitive 
to the value of c,, to make it clear that the earlier results are not com- 
patible with the present heat capacity measurements. The moduli 
given by McSkimin and Bond (1957) yield @(elastic)—2240+-5°K when 
de Launay’s (1956) tables are used for the calculation. This value 
agrees with ©) within the assigned uncertainties. It may be noted here 
that Burk and Friedberg (1957) give @,=2246°K, which was obtained 
by linear extrapolation of results for the region 45°<7'<100°K plotted 
in the form of @) against 7?. This would appear to correspond to a 
‘pseudo ’ a734-b7'> region (Barron and Morrison 1957) and hence would 
be expected to yield slightly too large a value of @). 

The good agreement between @, and @(elastic) suggests that McSkimin 
and Bond’s elastic constants are the most reliable, but it should be 
stressed that @ is a direct check on the value of c,, only ; the values 
of ¢,, and c,, need independent confirmation. In particular, new 
measurements of the compressibility of diamond would be valuable : 
the earlier measurements of Adams (1921) and of Williamson (1922) are 
inconsistent with the moduli given by McSkimin and Bond. 
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4.5. Comparison with Models 


Smith (1948) has computed an approximate frequency spectrum for 
diamond and has derived values of @,, at selected temperatures, Her 
results for the region 7'<300°K are plotted in fig. 1; it is evident that 
they depart seriously from @,, as calculated from ihe measured heat 
capacities. The force constants for the model were based in part upon 
the elastic constants of Bhagavantam and Bhimasenacher (1946), which 
now appear to be less reliable; this will account for some of the 
discrepancy. 

In addition, however, reservations may be held about the assumption 
of central forces for second neighbour interactions. It may prove useful 
to try in preference a five parameter theory assuming quite general forces 
between both first and second neighbours. This would neglect only 
interactions between pairs of atoms which have no common nearest 
neighbour. The force constants might be fixed by the elastic constants 
of McSkimin and Bond (1957), the first order Raman frequency shift 
(Robertson e¢ al. 1934) and ©,, the high temperature limiting value of 
@,. Unfortunately, existing heat capacity results at high temperatures 
(Magnus and Hodler 1926) do not seem to be sufficiently accurate to give 
a good estimate of 0,,. 


$5. SUMMARY 


The experimental results given here resolve some of the existing 
uncertainties about the thermal properties of diamond. In particular, 
they show that the form of the low frequency spectrum is not unlike 
that of other cubic crystals such as the alkali halides. The close agree- 
ment between @, and @(elastic) confirms the value of one of the elastic 
moduli (¢,,) measured by McSkimin and Bond (1957). 
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AP BND TX 


Measured Heat Capacities of Diamond 


(units—cal/g atm deg) 


TR hie Of TK) Op) PER Es ol! Ce) es 


12-833 | 0-000115 | 60-540 | 0-01035 | 120-283 | 0-1124 | 218-511 | 0-7146 
12-968 | 0-000126 | 62-189 | 0-01131 | 123-697 | 0-1242 | 222-107 | 0-7449 
16-015 | 0-000192 | 63-648 | 0-01237 | 126-975 | 0-1360 | 225-595 | 0-7761 
16-745 | 0-000190 | 65-392 | 0-01349 | 127-870 | 0-1385 | 229-264 | 0-8068 
18-631 | 0-000241 | 66-750 | 0-01433 | 131-154 | 0-1515 | 232-816 | 0-8404 
19-757 | 0-000327 | 68-591 | 0-01573 | 184-485 | 0-1654 | 236-364 | 0-8711 
21-304 | 0-000424 | 70-070 | 0-01681 | 137-805 | 0-1794 | 239-816 | 0-9033 
22-464 | 0-000501 | 71-787 | 0-01814 | 141-116 | 0-1944 | 243-176 | 0-9341 
24-100 | 0-000600 | 73-386 | 0-01965 | 144-495 | 0-2097 | 246-452 | 0-9644 
25-276 | 0-000708 | 74-995 | 0-02116 | 147-945 | 0-2271 | 249-159 | 0-9887 
26-994 | 0-000831 | 76-409 | 0-02220 | 151-444 | 0-2448 | 252-541 | 1-0220 
28-299 | 0-000975 | 76-685 | 0-02279 | 154-989 | 0-2641 | 256-039 | 1-0553 
29-992 | 0-00117 79-169 | 0-02546 | 156-186 | 0-2703 | 259-648 | 1-0926 
31-332 | 0-00135 79-452 | 0-02555 | 159-510 | 0-2886 | 263-168 | 1-1232 
33-407 | 0-00161 80-868 | 0-02762 | 162-867 | 0-3083 | 266-789 | 1-1614 
34-606 | 0-00178 81-972 | 0-02895 | 166-298 | 0-3273 | 270-507 | 1-1967 
37-316 | 0-0022 82-410 | 0-02954 | 169-777 | 0-3503 | 274-134 | 1-2332 
38-046 | 0-00243 84-808 | 0-03258 | 173-316 | 0-3728 | 277-675 | 1-2686 
41-319 | 0-00313 85-550 | 0:03364 | 176-730 | 0-3947 
41-325 | 0-00313 87-112 | 0:03574 | 180-042 | 0-4175 
44-491 | 0-00394 88-700 | 0-03838 | 182-415 | 0-4348 
45-228 | 0-00413 90-190 | 0-04058 | 185-742 | 0-4570 
47-570 | 0-00481 92-201 | 0-04360 | 189-157 | 0-4812 
48-985 | 0-00529 96-025 | 0-05077 | 192-641 | 0-5072 
50-450 | 0-00579 99-904 | 0-05812 | 196-182 | 0-5339 
52-045 | 0-00648 | 103-085 | 0-06462 | 199-611 | 0-5592 
53-980 | 0-00717 | 106-381 | 0-07311 | 203-157 | 0-5881 
55:597 | 0-00785 | 109-756 | 0-08149 | 206-809 | 0-6167 
57-333 | 0:00865 | 113-212 | 0-09073 | 210-341 | 0-6465 
58-957 | 0-00944 | 116-719 | 0-1009 | 215-014 | 0-6841 
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ABSTRACT 


Skyrme has recently discussed the use of a model in quantum mechanics. 
His method is applied to the case of the anharmonic vibrations of a crystal 
lattice, and compared with a previous treatment by the present author. Some 
remarks are added which give a new and more physical interpretation of the 
results of this earlier work. 


$ 1. INTRODUCTION 


THe dynamics of a perfect crystal lattice has nowadays little interest 
except when the forces are strongly non-linear, or the lattice supposed 
to be in interaction with some internal perturbing system. The former 
case was treated some years ago by the present author (Hooton 1955 a) 
on the basis of a ‘model’ suggested originally by Born (1951). 
Recently, Skyrme (1957) has discussed the use of model systems in 
quantum mechanics in a wider context, and the philosophy of his approach 
is exactly that used in the papers quoted above. Indeed, Skyrme first 
illustrates his method by the academic example of a one-dimensional 
anharmonic oscillator. The physically interesting case of the coupled 
anharmonic oscillations of a lattice was developed by Hooton (1955 a) 
from precisely the same model and it is instructive to compare the two 
methods. Some generalization of Skyrme’s formalism is of course 
required, since the lattice must be regarded as a system in thermal 
equilibrium. <A simple way of doing this is indicated below, although 
it is not suggested that this approach should be followed in practice. 
For the anharmonic lattice, which is a system of great complexity, the 
methods used by Born (1951) and Hooton (1955 a) offer much greater 
hope of a simple solution. Indeed, we shall give here an alternative 
and more direct physical argument which yields the previous results 
immediately. In the present context Skyrme’s terminology can be 
used to advantage to give a very clear statement of the method of 
approximation which the two approaches have, in fact, in common. 


§ 2. THE DynamicaL MopEL 


The method used by Hooton (1955 a) can be formulated briefly as 
follows. We introduce collective coordinates and momenta q, p (suffixes 
are omitted for simplicity) to describe the coupled vibrations of the lattice 
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particles; they are obtained by a linear canonical transformation of 
matrix (¢) from the displacements X—X°® and momenta P of the indi- 
vidual atoms. The centres of vibration X° are parameters which must 
be determined by a method which takes full account of the non-linearity 
of the system—the only reasonable choice is to define them as the mean 
positions (XJ, the lattice being treated as a system in thermal equili- 
brium at temperature 7. The transformation matrix (¢) will also be 
determined if we can find an approximate form for the coupled motion 
in the lattice. Were the forces linear, this coupled motion would of 
course be a superposition of independent harmonic modes. Born (1951) 
has suggested that as a model of the non-linear motion we might again 
consider a set of ‘normal modes’ but now with frequencies w (which 
appear as additional parameters) determined in a non-linear fashion. 
The chief problem is in fact to choose some characteristic of such a set 
of harmonic oscillations to impose on the actual motion, for which it 
must hold approximately if the model is to be at all realistic, thereby 
gaining equations which will determine the parameters introduced 
above. 
We write} the actual Hamiltonian : 


H=typt+Vig), <> 4 


where the potential U(X) = V(q) has arbitrary (polynomial) form. 
The model is written : 


h=tSp4+4ye% es (8) 


The philosophy that this method has in common with Skyrme is to 
suppose that the level scheme of h is the same as that of H save for an 
approximately constant displacement ¢, at least so far as the low levels 
are concerned. In fact in a lattice at normal temperature there will be 
comparatively little excitation above the ground state. For solid 
helium—to which the results of this method have already been applied 
(Hooton 1955 a)—we can deal practically with the zero-level alone ; 
« is then approximately 4H), the difference in energy as 77—0. It is 
also implied that the eigenfunctions %,,(q) and ¢,(q) of (1) and (2) are 
similar—again an approximation which obviously fails for highly excited 
states. Skyrme’s formal development depends on the approximate 
identification of these eigenfunctions, as we shall see below. It should 
be remembered that all these remarks can be expected to hold only when 
the parameters have been fixed so as to fit the model to the actual system, 
if this is indeed possible. 

The preceding philosophy can be expressed by saying that the operator 


~ 


H—h—e=V(q)—e,  V(qg)=V(q)—3> 0%, . . . (3) 


7 The notation used differs from that of Born (1951) and Hooton (1955 a) 
in an obvious fashion, and is chosen to agree with Skyrme to facilitate comparison 
with his paper, 
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is very small, having approximately degenerate eigenstates of eigenvalue 
zero. ‘This means that in forming any thermal average with the density 
matrix 
p(H)=exp (—H/kT)=exp (—«/kT)p(h) exp {—(V—«)/kT}. (4) 

we may approximately disregard the last factor—an approximation 
which effectively replaces the correct eigenfunctions y;,, with the model 
functions ¢,, in those states where this is most likely to be true. 

Moreover, the free energy (from which all thermal lattice properties 
follow) can be written approximately 


F=—KkT In irp(H)>—kT Inirp(h)t+e. . . . . (5) 
The eqn. (5) requires merely that the average value of the operator 


exp {—(V—e)/kT} should be unity ; to first order in the operator V—e« 
this at the same time determines e: 


e=(V),. Rae re 5 4) sft (6) 
(Here ¢f(7)) ,=trie(h)f(q) }/trp(h), the thermal average constructed accor- 
ding to the model Hamiltonian.) 
Thus : F2—kT In trp(h)+(V),, Be Seon earth) 


which was in fact the form derived by Hooton (1955 a) by the use of 
perturbation theory. Once the various parameters have been specified 
to fit the model to the actual system, (7) can be completely worked out. 


§ 3. Forma USE oF THE MODEL 


In Skyrme’s formulation, the model operator h is described in its own 
space of variables €, and a formal eigenvalue equation is derived for the 
operator H—h—e in the combined gq, € space. The (approximately 
degenerate) eigenfunctions, written F(q, €), are defined by the relation 


Bq, £)= Lana Don™(€), ee arate bee (| 
with arbitrary coefficients a,. These eigenfunctions are to be determined 
so as to express the degree of correspondence between the states %,, and 
¢,. The discussion of § 2 identifies g and €, which here implies writing 


F(q, £)=118(q—&) ; 7 Re ee ED) 
in view of (8), this means putting 4,—¢, (with a,—1) for all states. 
Before criticizing this procedure, it is worth noting that if (9) is inserted 
in Skyrme’s variational principle : 


af { F*@, Q(H—h— PF, )dqdé=0 . . (10) 
one gets formally for « the equation 
tr(V(q)—e)=0. Mer Geet. dL) 


This is just the eqn. (6) in the limit 7’— oo, since only the Boltzmann 
weighting factors are missing. Variation with respect to the other para- 
meters contained in the operator is not easy to interpret even formally, 
since derivatives of the 5-functions are involved. 

The fault in this use of Skyrme’s formulation is clear. Whereas 
Skyrme dealt with the static energy problem, we wish to deal with a 


DZ 
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thermodynamic problem. In § 2 a direct method of approximation was 
followed (in a single g-space) which used the thermal weighting factors of 
the model, exp {—¢,/k7'}, to cut off the excited states. In (9) and (11) 
all states have been given equal weight. Thus, if we wish to use Skyrme’s 
very elegant formulation in terms of the combined q, €-space, the variables 
must be kept distinct and some weighting introduced to limit the 
approximation ,<¢,, to those states where it is likely to succeed. 
Obviously for a system in thermal equilibrium one should introduce the 
Boltzmann factors in a manner analogous to that of the direct method— 
an additional assumption which extends Skyrme’s method to the more 
general problem. 

The clue to this lies in the degeneracy of (8): any eigenfunction of 
H—h should deliver us the same energy-difference «. We might put 
A,=8ym for each level %,, in turn; the function F’,,(¢, €)—¢n(V)bm*(E) is 
then itself an approximate eigenfunction of /7—h, and if used in (10) 
would yield an estimate of « (and the other parameters). Clearly this 
estimate is good only so long as the approximation %,,~¢,, is justified. 
To cut off the excited states we might determine the parameters by 
forming the equation 


(af [ Fat(g, \(H—h—F,,(q, €) dq dé),=0, . . (12) 


which to a first approximation minimizes the thermal average of the 
mean-square deviation of H from the model value h+-« with respect to 
the states ¢,. For «, (12) gives immediately the eqn. (6). 

The preceding extension of Skyrme’s arguments has been given to 
illustrate how the spirit of his method can be adapted’to deal with a 
system in thermal equilibrium, which poses problems of a different 
character to those arising in his own examples. However, it is not 
suggested that in the present case this approach would provide any 
useful advantage over that described in § 2. 


§ 4. DynamicaL DEFINITION OF THE LATTICR 


In the light of the preceding paragraphs, some remarks can usefully be 
added on the direct methods of fixing the parameters X°, (e) and w 
proposed in Born (1951) and Hooton (1955 a). 

We wish the mathematical centres of vibration X° to take account 
of the fact that the mean positions of the atoms may differ considerably 
from the configuration of minimum potential energy. As remarked 
above, the only physically reasonable assumption is to define the X° as 
the mean positions (X) themselves. The system is then to be thought 
of as vibrating about its mean configuration, and we have to determine 
the coordinates q in a way which best describes the similarity of this 
motion to the harmonic vibrations of the model. 

This definition of the X° was proposed by Born (1951) and is equivalent 
to writing 


(7) ORY weet A resent epee tray 
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it in fact imposes on the actual system (1) a trivial characteristic of the 
model (2), viz., that the coordinates have zero mean values. 

At the same time, the coordinates q were fully defined by imposing 
on the system a second characteristic of the model, viz., that different ¢ 
should be uncorrelated and should have the appropriate mean-square 


deviations : (94) =8yq'f(w, 7). Rieter eter ae te» 


Here f(@, T)= <q"), is a well-known thermal average for the harmonic 
oscillator. F 

According to a perturbation expansion of the density matrix (4) 
in power of q (Hooton, Edinburgh thesis) the eqns. (13), (14) lead piecemeal 
to the remarkably simple expressions : 

<dV/oq),=9, pin eee ome a8 (1D) 
{02V [0.0 9°) n= 8 49" re be A CG) 
The immediate significance of these equations will appear directly. 

In an attempt to get away from the small non-linearity implied in the 
expansion in powers of ¢ a variational method was proposed in Hooton 
(1955 a), the free energy being minimized with respect to the parameters 
introduced above. It should perhaps be emphasized that this is a useful 
method only if the variation is applied to the approximate form of F as 
given in (7). The approximation is then that expressed in (3) and (4), 
which treats the operator V(q)—e« as small, regardless of what powers of 
q it contains. The variational method then gives just the eqns. (15), 
(16) in full. But although it is thus compatible with (13) and (14) 
within the approximation (and (13) at least must be insisted upon), 
this approach does not bring out clearly which characteristics of the 
non-linear system it is proposed to identify with those of the model. 

A further point should be emphasized : all the preceding holds equally 
well for an arbitrary structure (say a single molecule) or a crystal lattice. 
The original intention of this work was to consider non-lattice solutions. 
This intention is no longer of interest, and we may take the mean con- 
figuration X° to be a regular lattice. In that case symmetry properties 
of the crystal should make the eqns. (15) redundant: the lattice (X°) is 
given, its parameters being determined ultimately from the free energy 
in the usual way. Lattice symmetry certainly makes the first term of 
each eqn. (15) identically zero; for the particular case of a one-dimen- 
sional chain the theory can be so developed that the next (cubic) term of 
(15) also vanishes identically (Hooton 1955 b). It seems reasonable to 
hope that this requirement can be satisfied in general, so that in effect 
only the eqns. (16) remain. It is, however, useful to treat both sets 
of eqns. (15), (16) together in what follows. 

A direct interpretation of (15), (16) can readily be found. If the 
system were harmonic with potential V’(q), equations of this type would 
express just the same characteristics of the normal modes as do (13) and 


ae (0V'/8q)=(8V'/8q),-o + o*(qy=wrX(qy=0  . . . (17) 
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(with the usual definition for the X°), and 
(02V'/0qaq' )=(02V"/0q09') 9, g’-0= © 0 aa" mag chee) 


(which diagonalizes the quadratic terms, thereby ensuring zero correla- 
tions). The vibrations of a harmonic system are in fact entirely character- 
ized by these equations, though the averages on the left are a rather 
formal way of writing things in the harmonic case. 

It is otherwise when the system is non-linear. We could well choose 
just these equations, which use the same fundamental characteristics 
of the harmonic model as are used in (13) and (14), and demand that they 
be (approximately) satisfied for the real vibrations : 


<dV/oq)=0, acy atts Ow syle! ae NLD 
(0?V/dqdq’ )=w 549°. > oti oe, goes 


We thus again define coordinates gq which should exhibit just those 
properties of the real vibrations that approximate to the characteristic 
properties of a set of normal modes. The method of approximation 
underlying the use of the model, as expressed in (3) and (4), may now 
be applied directly in evaluating (19) and (20) : we work out the averages 
approximately using the model Hamiltonian, so that (19), (20) go over 
into the previous eqns. (15), (16). 

Within the limits of the approximation we have achieved the same 
results as previously ; in particular, the approximate eqn. (19) is com- 
patible with (13), as we must insist if the reference configuration X° is 
to have a proper physical meaning. The advantage of this approach 
is to show very directly how the method depends on the use of the model. 
The single assumption is that the operator H—A—e can be made for all 
practical purposes very small. 
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ABSTRACT 


A three-fold counter telescope has been used to study the radiation 
penetrating 10 cm of lead at a latitude of 18° N. A barometric coefficient 
of —0-14+0-02% per millibar was deduced from the semi-diurnal variation. 
The average counting rate decreased on days of high geomagnetic activity 
and at the same time there was a large increase in the diurnal variation. 


§ 1. INTRODUCTION 


ALTHOUGH it is obviously desirable to obtain data on any geophysical 
problem from as many different parts of the world as possible there is an 
additional advantage to be gained from studying cosmic rays in the 
tropics. In temperate latitudes the variations of pressure, and other 
atmospheric variables, are large and irregular, whereas in the tropics the 
most noticeable feature of the pressure records is a semi-diurnal variation 
superimposed on a relatively steady background. There is therefore 
much less uncertainty in applying corrections to the observed intensities 
in order to determine variations in the primary cosmic ray intensity. 
The experiment described here was a relatively simple one but was 
sufficient to show the value of working at this site. 


§ 2. APPARATUS 


The three-fold counter telescope shown in fig. 1 had a coincidence rate 
of about 13000 per hour. Recordings were taken every half hour, 
although only the hourly totals were used for most of the analysis. The 
design of the equipment followed current practice in this field and parti- 
cular attention was paid to stabilizing the power supplies. Unfortunately 
the electricity supply was liable to interruptions caused by thunder- 
storms ; these were particularly prevalent during the summer months. 

This apparatus was operated under a light roof at an altitude of 180 
metres above sea level. Recordings were taken between September 1955 
and August 1956 but those for a number of days each month had to be 
discarded through electrical supply failures, apparatus faults or routine 
checking. 
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The geographic latitude of the recording site is 13 oe N and the 
longitude 76° 44’ W. Sea level meteorological data were obtained from 
the Palisadoes airport a few miles away. No upper air soundings are 
taken in Jamaica. Magnetic data were supplied by the San Juan 
Observatory in Puerto Rico. 


Fig. 1 
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Diagram of counter telescope. All measurements are in centimetres. 


§ 3. RESULTS 
3.1. Analysis of Daily Totals 


A regression analysis between the total number of counts each day and 
the daily average of atmospheric pressure gave a barometric coefficient of 
—0-228-+0-032%, per millibar. The statistical error here is large because 
the changes of pressure were always small ; the standard deviation of the 
daily averages was only 1-3 millibars. 

After correcting the daily totals for the pressure changes, using the 
above coefficient, the remaining variations were still much larger than 
the inherent statistical ones. These variations might partly be accounted 
for by other meteorological factors. Data for these were not available 
but, in view of the general regularity of tropical weather, it seems un- 
likely that their influence could be sufficient to explain the variations. 

Examination of the records suggested that the cosmic ray intensity was 
affected by the prevailing magnetic activity. This was first established 
by a partial regression analysis using the barometric pressure and the 
daily sums of the three hourly AK indices as the independent variables. 
The barometric coefficient obtained in this way had the same value as 
before and therefore it is permissible to correct the daily totals for the 
pressure variations and then plot the resulting daily averages against the 
K numbers. Figure 2 shows that the intensity decreases with increasing 
magnetic activity and that the dependence on the A number is not a 
linear one. 

3.2. Non-periodic Variations 

The 60°, increase in intensity during the solar flare of 23rd February 
1956 has already been reported (Barton and Stockhausen 1956). During 
the period for which observations were made it was not possible to dis- 
tinguish any ‘Forbush’ type decreases of intensity. In order to 


— 
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examine more carefully whether the magnetic phenomenon of sudden 
commencements influenced the cosmic ray intensity the times of these 
events were obtained from the San Juan Observatory and the cosmic ray 
readings superimposed for a number of events. Figure 3 (a) shows the 
values for two-hourly periods, after correcting for pressure variations, 
for a total of nineteen events. In eight cases, for which there was no 
second event interfering and for which the records were sufficiently 
complete, it was possible to study this behaviour over a longer period, 
as shown in fig. 3(b). It is apparent that, on the average, a sudden 
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Dependence of mean counting rate on magnetic activity. 


commencement leads to a small drop in intensity which continues for 
several days; it is regarded as fortuitous that the decrease is larger for 
the particular events selected for fig. 3(b). Since the days following 
sudden commencements are aften ones with high A numbers this result 
is essentially equivalent to that reported in § 3.1, although there is also 
an indication that the intensity may be slightly above its normal value 
immediately prior to a sudden commencement. 


3.3. Daily Variation 


The hourly totals for the number of cosmic rays were first grouped for 
each calendar month and then analysed for the diurnal and semi-diurnal 
components. The results for the different months were more variable 
than would be expected statistically but did not show any systematic 
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variation. A second analysis was therefore carried out in which the 
readings were grouped according to their K numbers. The results of 
this latter calculation are shown on two harmonic dials in fig. 4. Also 
plotted on these dials are the vectors for the daily variation of baro- 
metric pressure ; these are average values computed from the meteoro- 
logical records for the same period. 


Fig. 3 
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Averaged cosmic ray intensity at times of sudden commencements. 


Variations of the Cosmic Ray Intensity in Jamaica 59 
Fig. 4 
A DAYS WITH O<K< 10 


00h 


Banas « 10<K<20 
C hear 20<Ks 50 
De. » 30<K<40 


ERRORS MARKED ARE STANDARD DEVIATIONS 
OP IS PRESSURE VECTOR 


Gree 0.5 


0-4 0.5 


0.6 


18h 


I2h 


Ol 02 03%04 0.5 
09h 03h 


(0) 


Harmonic dials of pressure and cosmic ray intensity variations. 
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For days which are relatively quiet magnetically it is seen that the 
semi-diurnal variation is the more important one. For the average of all 
days with K <20 the amplitude of this harmonic is 0:15+-0-02 °/, and its 
maximum is at 4.00hr. It is almost opposite in phase to the second 
harmonic of the pressure variation and is presumably caused by this 
The necessary barometric coefficient would be —0-14+-0-02%, per millibar. 

The diurnal variation is very small for magnetically quiet days. For 
days with K<10 its phase is opposite to that of the pressure variation 
but the errors are too large for this result to be significant. All that can 
be said is that any residual diurnal variation, after correcting for pressure, 
must have an amplitude of less than 0-1%%. 

With increasing magnetic disturbance both diurnal and semi-diurnal 
components increase. This is particularly striking for the diurnal com- 
ponent which, after correcting for the pressure variation, has an amplitude 
of 0-60-40-08% on the most disturbed days. 


§ 4. Discussion 
4.1. Barometric Coefficient 

The value for this found from the semi-diurnal variation agrees closely 
with that found by Sarabhai, Desai and Kane (1953) from low altitude 
readings at geomagnetic latitude 13° N. The semi-diurnal pressure 
variation is generally believed to be explicable in terms of an atmospheric 
solar tide. According to the calculation of Weekes and Wilkes (1947) the 
relative amplitude and phase of the pressure variation are substantially 
constant from sea level up to a height of 25 km. Since this covers most 
of the region which is important for the production and decay of mesons 
the theory of Janossy (1948) for the pressure coefficient should be adequate. 
Putting in values obtained from Rossi’s (1948) curve for the absorption 
of the hard component near the equator this theory predicts a coefficient 
of —0-17°% per millibar. 

The value obtained from the day to day changes of counting rate 
appears to be slightly larger than that deduced from the semi-diurnal 
variation. Although this is not yet established with certainty it is 
possible that temperature effects may lead to a real difference. Com- 
parable differences between coefficients determined by using daily and 
monthly averages have been reported by Dolbear and Elliot (1951). 


4.2. Geomagnetic Effects 

The harmonic dials of fig. 4 show that the degree of magnetic activity 
plays a dominating role in the daily variation. On magnetically quiet 
days the results can be interpreted directly in terms of radiation which 
is incident isotropically, within 0-1°,, and modulated solely by atmo- 
spheric effects. With increasing magnetic activity the daily variation 
increases very markedly. The diurnal component exhibits the more 
striking increase and at the same time shifts in phase to a value almost 
opposite to that prevailing on quiet days. (The results plotted in fig. 4 
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are average ones; there is some indication that individual days show 
appreciably different behaviour.) Combining this result with the de- 
crease of average counting rate on disturbed days, which has about the 
same magnitude as the amplitude of the diurnal variation, it seems 
natural to conclude that the physical process is one in which a proportion 
of the cosmic ray beam is removed. This effect reaches its maximum in 
the early evening. It is not possible to state definitely that the cosmic 
ray intensity never exceeds its mean value at any time throughout the 
day but the results certainly show that reductions in the intensity are 
predominant. 

The great influence of magnetic activity on the form of the daily 
variation was pointed out by Sekido and Yoshida (1950) and Elliot and 
Dolbear (1951). The former authors found that the amplitude of the 
diurnal component at stations in Japan was 0-35°% during intense 
magnetic storms. Sekido and Kodama (1952) also noted that the 
phase of this component advanced as it became larger, in agreement with 
the behaviour shown in fig. 4. 

Unfortunately most analyses of the daily variation have been made 
without distinguishing between days of different magnetic activity. 
However, Ehmert and Sittkus (1951) have given results for days influenced 
by strong magnetic storms ; for a counter apparatus operating at 48° N 
they found that there was a reduction in the intensity during the latter 
part of the day. The magnitude of this reduction was 0:74°%, which 
would be equivalent to a diurnal component of amplitude not greater than 
half of this figure. Firor et al. (1954), using neutron counters at high 
altitude, found a peak-to-peak amplitude of 1-0°%, at the equator and 
1-4°%, at 48° N for days with a mean K value of 33:8. Comparison between 
results obtained under different conditions may be misleading but it does 
appear that the magnitude of the additional daily variation is not greatly 
dependent on the mean energy of the primary cosmic rays involved. It 
is also important to note that the minimum intensity is always recorded 
in the evening when using local time. 

Failure to separate days of different magnetic activity may well account 
for some of the anomalous results which have been reported. Thus 
Possener and van Heerden (1956) reported a reversal of the phase of the 
diurnal component during certain months. Although they studied the 
intensity at a time of minimum solar activity it is quite possible that 
there were sufficient disturbed days to mask completely the true character 
of the variation for undisturbed days. Similarly the change of the 
diurnal variation which Thambyahpillai and Elliot (1953), and Sarabhai 
and Kane (1953) found to be associated with the solar cycle may merely 
be a consequence of the varying proportion of magnetically disturbed days. 

The explanation of the enhanced daily variation is not known. Firor 
et al. have shown that it is very difficult to explain the results by means 
of the direct effect of the earth’s field on the cosmic radiation. A similar 
conclusion has been reached from studies of the 27-day variation by 
Simpson (1954) and van Heerden and Thambyahpillai (1955). On the 
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other hand, Elliot and Rothwell (1956), by analysing the variation observed 
in inclined directions, concluded that the variation is due to modulation 
of isotropically incident primaries by fields close to the earth rather than 
to an anisotropic primary distribution. These views are not necessarily at 
variance since the arguments of Firor et al. do not exclude an effect due 
to the earth’s field of about 0-1°%. Just as the behaviour of the earth’s 
field on quiet days is qualitatively different from its behaviour on dis- 
turbed days so it is possible that variations in the earth’s field are sufficient 
to explain Elliot and Rothwell’s results but are less significant for the 
effects which predominate on disturbed days. The assumption made 
here is that on the disturbed days both the earth’s field and the distri- 
bution of cosmic rays are influenced by some third factor, such as 
turbulent magnetic fields ejected from the sun and extending over much 
of the solar system. There is little direct evidence for such a picture 
but at least it is not contradicted by the results yet available. 


§ 5. CONCLUSIONS 


1. The barometer coefficient for the hard component of the cosmic 
radiation in Jamaica is —0-14+0-02°%, when determined from the semi- 
diurnal effect and —0-23--0-03%, when determined from day to day 
variations. 


2. On magnetically quiet days the diurnal and semi-diurnal compo- 
nents of the daily variation can be explained solely in terms of pressure 
changes. 


3. On magnetically disturbed days the mean intensity falls and there 
is a great increase in the daily variation. 
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A SINGLE crystal of zinc having a predetermined orientation was prepared 
from 99-999+°%, grade material by the Bridgeman technique and sub- 
sequently etched in a solution containing 200 g CrO., 15 g Na,SO,, 50 cm? 
of concentrated HNO, and water to make 1000 em? (Vinaver and Dreulle 
1955). 

On a face nearly parallel to the basal plane, several complex spirals 
were observed after etching for one hour. A typical example is given in 
fig. 1 (a), Pl. 1. These spirals were both left- and right-handed, and 
were surrounded by more finely etched areas. 

The pattern of fig. 1 (6), Pl. 1 was observed only in one location, in an 
area which was otherwise free from any ‘structure’. Interferometric 
observations revealed that the background was levelled, and that the 
‘loops’ stood about 0-15 micron above this background. The typical 
pattern is similar in aspect to an operative Frank—Read source (Frank 
and Read 1950). However, some of the detailed features, namely, the 
two ‘ dots’ near the centre of the pattern and the central ‘dip’ of the 
innermost line, are unexpected. 

An explanation for this peculiar and heretofore unreported etching 
behaviour is suggested, assuming that regions around the dislocation 
lines are saturated with respect to a solute element much more noble than 
zine, present at a very low overall concentration. A galvanic action 
would, therefore, occur, the material surrounding the dislocation lines 
dissolving at a lower rate than the rest. It should be emphasized, more- 
over, that the etching reagent was always found very discriminative with 
respect to structural features when the solution was not agitated, where- 
as it appeared to be an excellent levelling and polishing reagent when 
continuously agitated. 

Interpreting the pattern as a Frank—Read source, the length L can be 
measured (2 x 10° A), and the stress required to place the source in operation 
can be computed, based on Mott’s formula (Mott 1952) 

stress=Gb/L. 

Assuming 4:9 x 10° p.s.i. for G (Smithells 1955) and of 2-659 4 for b, the 
stress appears to be about 65 p.s.i., which is a reasonable value. 
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ABSTRACT 


The rolling texture of drawn polythene is interpreted in terms of slip and 
twinning processes characteristic of ordinary crystal plasticity in combination 
with the rubbery elasticity characteristic of macromolecular amorphous 
substances. The principal crystal plasticity elements operative are (100) [010] 
slip and (310) [130] and/or (110) [110] twinning. The observations appear 
to indicate that (110) [110] slip is a less easy mode of deformation than these, 
at least initially. Under more general deformation conditions, it appears 
that these modes predominate at small strain, tending to a texture with 
[010] aligned to the axis of extension, giving place to [001] slip processes, which 
align [001] and the molecular chains along the axis of extension when the 
strain is large. Under certain conditions a second crystalline phase also 
appears which in terms of general crystal plasticity would correspond to 
a martensitic type transformation and might arise through a regular sequence 
of stacking faults on {110} planes. 


§ 1. EXPERIMENTAL 


THe material of our experiments was melt pressed polythene film, 
Alkathene grade 2, about 0-1-0-15mm thick. Strips of this film were 
passed between rollers, the spacing between the rollers being adjusted as 
required. In another type of preparation strips of the film were first 
drawn by hand past the stage of necking to an elongation of 450%. 
Subsequently these drawn strips were passed between rollers, the direction 
of advance being along the length of the strip. Relaxation experiments 
were carried out by placing the drawn, or drawn and rolled strips on a 
temperature gradient bar at the required temperature. 

When a strip was mildly rolled the first sign of orientation was of the 
type where (100) was parallel to the plane of the film and [001] was parallel 
to the direction of advance, the orientation itself being only poorly defined. 
On increasing rolling this orientation became better defined but in addition 
to it a second type of orientation with (110) parallel to the film also 
appeared (fig. 1, Pl. 2). In the sample of fig. 1, Pl. 2, the length of the 
film increased by 11%, the width by 4% ; and thicknesses decreased by 
14%. The second type of orientation with (110) parallel to the film 
became more pronounced with stronger rolling, and predominantly this 
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orientation was produced when the drawn specimen was rolled (fig. 2, 
Pl. 2). (The effect of drawing was to orient [001] first.) The dimensional 
changes in the case of fig. 2, Pl. 2, were: increase in length 13-5%, 
increase in width 12-7°%, reduction in thickness 25°%. On mild rolling of 
the drawn samples there were faint signs of the (100) plane orientation 
besides the strong (110) orientation. Summing up, these results indicate 
that mild rolling produces preferentially (100), strong rolling (110) orien- 
tation, the relative proportion of the two orientations varying according 
to whether the sample was previously drawn or not. These observations 
agree essentially with those of Point (1953). 

Drawn and rolled specimens which were only slightly relaxed by an 
anneal at ~100°c changed from the original (110) to (100) plane orien- 
tation, [001] staying parallel to the direction of original draw (fig. 3, Pl. 3). 
The effect of stronger annealing will not be discussed here. It will only 
be recalled that on annealing of a drawn sample [001] gradually tilts away 
from the original draw direction while the position of [010] stays un- 
altered, until [100] becomes parallel to the direction of original draw (e.g. 
Keller 1955, Pl. Ic). Redrawing of such a sample has been examined by 
one of us earlier (Keller 1955) and it was reported that the crystals first 
turn around [001] prior to [001] alignment. On re-examining this effect 
we now find that this rotation around [001] is not continuous but results 
immediately in the orientation shown by fig. 4, Pl. 3. 

A further feature which became apparent is a new reflection at 4:59 A 
(also reported by others in somewhat different connection ; see later). 
As seen this reflection is present in figs. 1 and 2, Pl. 2, but is absent in 
figs. 3 and 4, Pl. 3. Its intensity relative to the others varied from 
sample to sample. 


§ 2. Discussion 


2.1. Introductory 


‘ Drawing’ chain-macromolecular substances, i.e. subjecting them to 
permanent extension of the order of unity or greater, whether they are 
crystalline or amorphous, produces a texture with preferred orientation 
of the molecular chains along the axis of extension. Crystalline polymers 
show different and rather complex orientation phenomena for smaller 
strains : we defer consideration of these till later in the paper. If the pre- 
drawn crystalline polymer is subjected to transverse compression, in- 
ducing only transverse displacement of the molecules, one may anticipate 
that the behaviour should be more nearly similar to that of non-macro- 
molecular substances, governed by familiar principles of crystal plasticity : 
except that, since in any attainable condition of the majority of polymers, 
the chains are neither completely aligned nor fully extended, some 
restoring stresses of rubbery nature will be present. Since rubbery 
elastic moduli are relatively weak, these latter stresses should have little 
effect until the strain is fairly large. 


P.M, 1} 


66 F.C. Frank, A. Keller and A. O’Connor on the 


In the light of these considerations we approach the interpretation of 
the rolling texture of drawn polythene, regarding rolling as sufficiently 
nearly the same as transverse compression. In those of our experiments 
where the polymer was not pre-drawn to begin with, conditions will not 
be very different. Here drawing took place during rolling, consequently 
at some stage compression will act on the material which is already drawn. 
The crystal structure (Bunn 1939) corresponds to the orthorhombic form 
of long chain paraffins (Miller 1928) illustrated in figs. 5 and 6. 


Fig. 5 

Yoo), _§_ pf 
K 
as 


(110) Twin. (170) Twin. 


a 


(310) Twin. — (370)Twin. 


Illustration of the effect of twinning on the orientation of the unit cell when 
viewed in the ¢ projection. 


2.2. Hapected Modes of Deformation 

The expected modes of crystal plasticity involving displacements 
transverse to the chain axes may be simply derived by inspection of the 
crystal structure. The shortest lattice translation (other than that along 
the c-axis) is 6 [010]. We therefore expect a slip-mode involving dis- 
locations with this as Burgers vector. Gliding dislocations should not 
be able to cut across the covalently bonded chain axes, so that the only 
possible slip plane for this dislocation is (100). The next largest trans- 
lation is a [100], with (010) as its necessary slip plane. As it is a con- 
siderably longer displacement, we should expect this dislocation to be 
less mobile. Since it is orthogonal to the [010] (100) slip system, it is 
always subject to the same driving stresses, and we should generally 
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expect the [010] (100) slip to relieve these stresses before | 100] dislocations 
move. [a6 0] dislocations would be expected to dissociate into pairs of 
partial dislocations with Burgers vector [}a, $b, y] and [4a, 4b, —y], y 
being an unknown quantity, each having a slightly shorter Burgers vector 
than 6 [010] if y is zero. These pairs would glide on (110). This is a 
reasonable slip system, but may not be observed because it is in com- 
petition with the twinning modes. 

Expected twinning modes are most simply deduced from Mallard’s law 
(cf. Frank 1953). The structure is orthorhombic pseudohexagonal. 
Hexagonal symmetry planes lost in the reduction of symmetry to ortho- 
rhombic are, in orthorhombic indices, (110), (110), (310) and (310). We 
anticipate these as twinning modes of the first kind: i.e. such that twin 
and matrix are reflections of each other in these planes as interfaces. 
These four modes make two conjugate pairs (310) with (110), (310) with 
(110), being pairs of planes orthogonal to each other in the related structure 
of higher symmetry. Hence no further modes are predicted from these 
considerations. In full, the corresponding twinning elements are : 


Composition Conjugate Shear Conjugate Semi- Shear 
plane plane direction axis angle S= 
K, Ky 1 Ne d 2 cot 2 
(110) (310) [110] [130] 41°27’ 0-249 
(310) (110) [130] [110] 41° 27’ 0-249 


The values of ¢ and S are calculated from the lattice parameters of Bunn 
(1939), a=7:40, b=4-96. 

Figure 5 shows the unit cell in ¢ projection when twinned according 
to these two twin laws. Figure 6 shows the structure in the case of (310) 
twinning. As seen a rearrangement of the planes of the carbon zigzags 
is also required. 

These expectations are in accordance with the observations of Koolvoort 
(1938) who observed (110) and (310) twinning in the normal paraffin 
C,, H;9. both as a result of deformation in the orthorhombic modification, 
and as a result of transformation to this form from the hexagonal modi- 
fication which exists near the melting point. Incidentally, his observa- 
tions are remarkably similar to those which may be made on corresponding 
metal systems. 

The appropriate stresses to operate any of these modes are rather 
similar, namely compressive stresses in a range of directions near to 
[010] or tensile stresses in a range near to [100], with a preference for one 
or other of the conjugate pairs according as the axis of principal stress 
deviates to right or left from these directions. Whether (310) or (110) 
twinning will operate more easily cannot be certainly decided on a priori 
evidence, but on weak empirical grounds (cf. Frank 1953) one could expect 
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the system of higher indices, i.e. (310), to be preferred. This appears to 
be supported by the present observations, but the (110) twinning may 
occur as well: since it produces about the same lattice orientation, 
it is neither easy to distinguish by x-ray diffraction nor of great importance 
for the plastic behaviour, whether both or only one of these modes occurs. 


Fig. 6 


(310) 


K, (310) 


(310) twinning in terms of the structure, in ¢ projection. The full lines show 
the parent, the dotted lines the twinned lattice. A, is the twinning 
plane common to both, K, and K,’ are the conjugate planes, which 
remain unaltered during twinning, in the parent and the twinned lattice 
respectively. Arrows show displacements of molecules in the twinnin 
shear. It is assumed that molecules lying on the composition plane EK, 
take up the symmetrical position parallel and perpendicular to the lane 
indicated in the drawing. 4 


») € ° . 
2.3. Consequent Deformation Behaviour 


On this basis we expect the following behaviour : 

Crystallites with (100) and (010) planes near to 45° to the axis of 
compression will deform by (100) [010] slip, undergoing a rotation which 
brings (100) more nearly normal to the compression axis. 

Crystallites with (010) planes near to being normal to the compression axis 
will twin, say on (310). This makes a moderate deformation (25 % cheaet 
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accompanied by a large change of orientation, by which the planes (100), 
(010), (110), (110), (310) and (310) undergo small rotations (zero for (310), 
of course) and transform into (110), (310), (100), (110), (310) and (010) 
respectively. 

The transformations of planes are the same for (110) twinning: only 
the associated small rotations are different. This is now a suitable 
orientation for further slip on (100) bringing that plane nearer to the plane 
of compression. 

When this orientation is closely approached, (100) slip should cease. 
No twinning is possible. The only likely mode of deformation in slip on 
(110). If the flow stress for this mode is larger than that for (100) slip, 
there will be no large departure from (100) orientation before (100) slip 
starts again. 

The outcome. for crystallites of any initial orientation, is a two 
component texture with (100) poles on either side of the pole or the rolling 
plane. If the flow stress for (110) slip is large, these (100) poles will be 
close together. In any case, allowing for the broad ares of indefinite 
orientation usually found with polymers, including ordinary polythenes, 
this would appear as a single arced (100) texture. This agrees with 
observation. 

We have remarked however that with larger amounts of deformation 
we must expect consequences from the fact that we are dealing with 
long chain molecules, imperfectly aligned and not fully extended. 
This refers both to molecules which may lie wholly or partly in an 
amorphous phase, and to molecules which may lie in folded configurations 
in the crystalline phase. The presence of such molecules must generate 
restoring stresses of a broadly rubber-like nature : that is, of low modulus, 
but persisting and increasing up to very large strains. These rubber-like 
stresses remain after the roller-pressure is removed, being equivalent to 
sideway compression or traction normal to the rolling plane. Their 
effect on the oriented crystallites should be to induce twinning, whereby 
the (100) planes, parallel to the rolling plane, rotate slightly and trans- 
form into (110) planes. This will produce exactly the two-component 
texture, both components having (110) planes preferentially parallel to 
the rolling plane, that is observed. This can be seen directly from fig. 5 
if the plane of rolling is taken to be parallel to the (100) planes in the 
untwinned crystal. 

In this texture, in which various crystallites, or parts of a crystallite, 
have twinned in alternative ways, elastically stored energy persists 
around the edges of every twin, balanced against a reduction in the energy 
associated with the rubber-like stresses. With a gentle anneal, insufficient 
to melt the crystals, the latter stresses can relax by molecular migration, 
as is regularly observed for rubber-like stresses in a non-crosslinked 
polymer. The crystallites can then assume a state of lower energy by 
twinning back again, so that the texture reverts to the single, or nearly 
single, (100) texture, as observed. 
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According to the foregoing, compression along [010] would produce 
twinning. Consequently the same is expected for tension along the | 
perpendicular direction [100]. This in fact is observed in the case of 
specimens which were drawn, relaxed and subsequently redrawn (fig. 4). 

The two possible modes of twinning can only be distinguished by small 
differences in the x-ray pattern. The ares being broad it is difficult to 
locate them with sufficient accuracy. However it was possible to reach 
a conclusion by comparing the angular positions of the 200 and 110 
along their respective circles. It can be deduced (see fig. 5) that the 200 
reflection should be nearer to the equator in the case of (110) than in that 
of (310) twinning. Accordingly we appear to have (310) twinning in fig. 2 
and (110) twinning in fig. 4. These modes of twinning were found to be 
characteristic of the two types of samples. When reaching these con- 
clusions it was assumed that before twinning the texture was closely 
single (100). If the initial texture was imperfectly aligned, then crystal- 
lites in orientations deviating in opposite directions would exhibit pref- 
erence in twinning between (110) and (110) modes, or between (310) and 
(310) modes, in either case displacing the (200) poles further from the 
meridian. There is evidence of such displacements by selection, the 
observed angles from meridian to the mid-points of (200) ares being 59° 
for the first type of specimen compared with 534° calculated for (310) 
twinning from the ideal orientation, and 68° for the second type of speci- 
men compared with 664° calculated for (110) twinning from the ideal 
orientation. 

The assumption, involved in the foregoing interpretation, of a compar- 
atively high critical stress for slip on {110}, is not one which could be 
made with confidence a priori. The energy of the stacking fault between 
two partial dislocations for sip on {110} is an unknown quantity. If 
it is sufficiently low, the critical stress for this slip mode could be about the 
same as, or even less than that for slip on (100) ; and no theoretical basis 
exists for predicting the relative magnitude of critical stresses for slip and 
twinning processes. However, the alternative assumption of relatively 
easy (110} slip does not appear capable of explaining the observations. 
If this process occurred in preference to one of the twinning processes 
under the stress which favours these, i.e. compression along an axis in 
the neighbourhood of [010], it would produce a strong component in the 
texture having (010) normal to the compression axis. This would be a 
stable orientation by reason of alternate slip on (110) and (110), the 
compression axis lying in the obtuse-angled sector between these two planes 
This is in disagreement with the observations. We infer that twinning 
does occur in preference to {110} slip. 

For a compression axis near to [100] the twinning processes do not 
enter into consideration. The (100) orientation (which is observed 
to occur as a texture component) would be unstable, the compression 
axis now lying in the acute-angled sector between (110) and (110). There 
should be two other stable orientations, one lying between (100) and (110), 
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the other between (100) and (110), in the obtuse sector in each case. We 
assumed above that (100) slip was easier than {110} slip, so that these 
two stable orientations are close to (100) and can be identified together 
as the observed (100) component of texture. Assuming, alternatively, 
that {110} slip is much easier, these two stable orientations would be close 
to (110) and (110). This could be identified with the {110} texture 
observed after strong compression but we should be left without an inter- 
pretation of the reversion to (100) texture on annealing or for the (100) 
component appearing after light rolling. 

A more complex possibility which cannot be rejected on the evidence 
is that, while (100) slip has a lower critical stress than {110} slip, it has a 
higher rate of work-hardening, so that the (100) texture develops first, 
and then changes to the {110} texture when this becomes the easier slip 
mode. The reversion to (100) texture on annealing remains unexplained 
in this case but might be attributed to an oriented recrystallization. 


§ 3. Exrra REFLECTIONS AND STRESS INDUCED PHASE TRANSFORMATION 


In the course of our work on polyethylene we have noticed reflections 
which cannot be indexed on the basis of the orthorhombic unit cell. 
In the rolled or pressed specimens in particular a strong reflection appeared 
corresponding to a spacing of 4:59 a, the corresponding lattice planes 
being parallel to the plane of rolling or pressing. The relative intensity 
of this reflection was variable which indicates the presence of a second 
crystalline phase. We also found that associated with this there are other 
reflections with correspondingly varying intensity which are less evident 
because of blending with other principal reflections. Such observations 
are not new. The reflection at about 4:5 A has been reported previously 
(e.g. Pierce et al. 1952) and a number of other new reflections have also 
been frequently quoted. Quite recently a new unit cell has been suggested 
by Teare and Holmes (1957) which would be a triclinic polymorph of the 
usual orthorhombic cell in analogy with similar polymorphism in paraffins 
(e.g. V. Sydow 1956). We concur with the principles of the interpretation 
of Teare and Holmes, but on the basis of our data we cannot accept their 
proposed unit cell. Our findings will be published separately once the 
work still in progress is complete. At this place only one aspect of 
this polymorphism will be discussed. 

It is a common feature of all observations that the strong reflection 
at about 4-5 A appears in specimens which were severely deformed. It 
appears, therefore, that we are probably dealing with a stress-induced 
phase transformation to an otherwise unstable crystal modification. Two 
important distinctions must be drawn in classifying stress-induced crystal 
transformations. Firstly, the stress may change the relative thermo- 
dynamic stabilities of two phases, or merely provide the ‘ activation ’ for 
a transformation to a stabler phase. Secondly, the mechanism of 
transformation may be incoherent, in which case hydrostatic pressure 
can be the only part of the stress relevant to the change in relative stability, 
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and change in volume the only thermodynamically relevant characteristic 
of the transformation itself; or alternatively, it may be coherent, pro- 
ceeding like a mechanical twinning process, in which case shear stresses 
will also contribute to the change in relative stability. Coherent * shear 
transformations’ are familiar in metals, and are called ‘ martensitic ’ 
after the classical example occurring in steel. The polymorphism of the 
normal paraffins and their derivatives, in crystal structures which essenti- 
ally differ only by shears and by molecular rotations about the chain 
axes, suggests that they should be particularly apt for martensitic trans- 
formation. The thermodynamic requirement for such a transformation 
to a normally unstable phase is that the product of strain in transformation 
and corresponding stress should exceed the difference in free energies per 
unit volume of the two unstressed phases. Since, in polymers, the 
difference in free energy between alternative crystal structures is 
governed only by van der Waals forces, whereas covalent bonds contribute 
to the ultimate strength, which limits the applicable shear stresses, there 
is a relatively high a priori likelihood that this condition may exist, 
compared with the case of metals in which the same kind of interatomic 
forces govern both quantities. The geometric requirement is that one of 
the three principal strains in the transformation be zero, while the other 
two are of opposite sign. If there is no volume change, this means simply 
that the deformation is a shear. There then exist necessarily two 
planes, both in general irrational, which are subject to no deformation, 
and either of which may serve as the planar interface across which the 
transformation takes place. Qualitatively, such a transformation closely 
resembles mechanical twinning. 

A particular hypothetical mechanism for a stress-induced phase- 
transformation would be the dissociation of pairs of partial dislocations 
on the {110} planes. In the first instance this produces single stacking 
faults on these planes. With a screw dislocation intersecting the plane 
it may produce a regular sequence of stacking faults constituting either a 
twinned crystal or a new phase. Alternatively, however, an irregular 
sequence of stacking faults of sufficiently high density might be respon- 
sible for the anomalous reflections. Unfortunately, though it is a very 
reasonable presumption that stacking faults of this kind can exist, the 
available degrees of freedom in molecular rotation about, and translation 
along, the c-axis prevent any reliable presumption of their detailed 
structure. All we are able to say is that a mechanism of this kind is 
reasonably likely, and would cause this sort of transformation to accom- 
pany extensive slip on {110} planes. 


§4. DEFORMATION UNDER MORE GENERAL CONDITIONS 
General empirical experience regarding crystal plasticity leads one to 
expect easiest slip along the directions of shortest lattice translation. 
This is reasonable, because the corresponding dislocations have the 
smallest Burgers vectors. The square of the Burgers vector enters into 
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the expression for dislocation energy, and energy barriers for dislocation 
motion are expected to be broadly proportional to their self-energy. By 
this rule, the easiest modes of slip in polythene should have the slip 
direction [001], for which the lattice translation ¢c=2-5 A, compared with 
a=74A, At first sight, this agrees with the broad experience that 
‘drawing’ polythene aligns the c-axes to the draw-axis. We ought, 
however, to be cautious, because we are dealing with elastically very aniso- 
tropic crystals. Stretching along the c-axis is opposed by covalent bonds : 
transverse stretching only by Van der Waals forces. This must raise the 
energy of the dislocations with Burgers vectors not in the a 6 plane. 
Compression along the c-axis, of the magnitude associated with dis- 
locations, is liable to be taken up by chain folding, which may allow the 
energy to be lower, but is likely to increase the energy barrier for motion, 
besides making it difficult to estimate. Thus we have the alternative 
theoretical possibility that slip along the b-axis, [010], may be easiest. 
The consequence would be that for moderate strains the b-axis would 
approach the principal tensile axis. This would make an orientation 
stable against twinning under that stress, and crystallites initially far 
from that orientation could make a substantial discontinuous approach 
to it by twinning. A fibre texture with axial b, and transverse a would 
develop. The deforming crystallites would also tend to acquire trans- 
verse c-orientation, but crystallites with near-axial c-orientation will not 
deform except by secondary processes, so that a less complete preferred 
orientation of c-poles should be produced. At a later stage, the chain 
continuity should enforce a new preferred orientation with axial c. A 
variety of observations are in agreement with this picture. For example : 
4.1. Drawing Behaviour 

Drawing produces [001] orientation ; the c-axes, the direction of the 
molecular chains, align in the draw direction. This chain alignment, 
however, is not gradual. In the initial stages of slow drawing at first the 
[010] poles congregate in belts more or less near to the draw axis (e.g. 
Keller 1955). The exact position of this accumulation of [010] poles 
varies with the conditions of drawing and with the type of specimen. 
The alignment of [001] starts only after this stage. In any interpretation 
of the drawing behaviour the spherulitic texture of the material would 
have to be considered (Keller 1955). This is beyond the scope of the 
present paper. It is apparent, however, that the presence of a [010] slip 
will contribute to the initial tendency of the b axes to move towards the 
draw axis. 

4.2. Relaxation Behaviour 

In highly drawn materials, [001] poles are approximately axial, [100] 
and [010] poles approximately equatorial. Heating produces a large 
longitudinal contraction, during which [010] poles retain their equatorial 
position, while [100] poles approach the axial, and [001] poles the equa- 
torial position (e.g. Keller 1955). We do not completely understand this 
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process. It is simple enough if we regard the material as crystallites 
embedded in a rubbery matrix, without intertransformation between 
crystalline and amorphous phases. We do not know to what extent the 
latter occurs, nor, if it does not occur, what prevents it. If it occurred freely, 
we should expect the material to crystallize in the drawn configuration, 
with axial [001] and even an elongation instead of a contraction. How- 
ever, accepting the fact of contraction, equivalent to axial compression, 
and supposing that crystallites largely preserve their identity during the 
relaxation, the stability of the equatorial position of (010) poles will follow 
from the supposed ease of [010] slip, and the orientation to which crystal- 
lites should tend is that with the slip plane (100) normal to the axis of 
compression. This is in accordance with the observations. 


4.3. Redrawing of Relaxed Specomens 
If the material relaxed to the state in which [100] is substantially 
axial is redrawn, (100) poles are displaced from the meridional position 
by about 67° without assuming intermediate positions (fig. 4). As 
discussed above this is what we should expect, since all crystallites are 
favourably oriented for twinning. [001] alignment sets in progressively 
after this stage. 


4.4. Orrentation within the Spherulites 


The crystallites on any radius of aspherulite have [010] parallel to this 
radius. This suggests what may be a profitable approach to the unsolved 
problem of how the spherulitic texture in polythene (and other polymers) 
originates. The foregoing considerations suggest that we may generally 
regard an axis of [010] textural orientation as indicative of the direction 
of a moderate extension. If this is applicable to the genesis of the 
spherulites it implies that crystallization is followed by contraction in 
towards a centre. It is not suggested that this will account in full for 
morphology of the spherulites. 


REFERENCES 
Bunn, C. W., 1939, Trans. Faraday Soc., 35, 482. 
Frank, F. C., 1953, Acta Metaill., 1, 71. 
Keer, A., 1955, J. Polym. Sci., 15, 31. 
Koxvoort, EH. C. H., 1938, J. Instn Petrol. Tech., 24, 338. 
Muximr, A., 1938, Proc. roy. Soc. A, 120, 437. 
Prercz, R. H., Torpena, J. P., and Bryant, W. M. D., 1952, J. Amer. chem. 
Soc., T4, 282. 
Port, J. J., 1958, J. Chim. phys., p. 76. 
Sypow, E. v., 1956, Ark. Kemi, 9, 231. 
TEaRE, P. W., and Hormus, D. R., 1957, J. Polym. Sci., 24, 496. 


Prom 


Charged Dislocations and the Strength of Ionic Crystals} 


By J. D. Esuensy, C. W. A. Newey and P. L. Prarr 
Department of Physical Metallurgy, University of Birmingham 


and A. B. Lipiarp 
Department of Physics, University of Reading 


[Received November 13, 1957| 


ABSTRACT 

If the energies required to form positive and negative ion vacancies in 
an ionic crystal are unequal, then in thermal equilibrium dislocations in the 
crystal will be electrically charged and surrounded by a Debye—Hiickel cloud 
of vacancies. If the vacancy cloud is immobile a finite force is required to 
separate the dislocation from the cloud, and so the crystal will possess a 
yield stress below which plastic flow will not occur. The presence of divalent 
unpurities modifies the magnitude of the charge on a dislocation, and may 
even reverse it. If precipitation of the impurity or association of impurity 
atoms and vacancies can occur, the concentration of impurities may be a 
complicated function of temperature. The yield stress of the crystal may 
then exhibit maxima and minima when plotted as a function of temperature. 
Hixperimental results showing this behaviour are presented and tentatively 
compared with the theory. 


$1. [NTRODUCTION 


THE number of Schottky defects (vacant anion and cation sites) in an 
ionic crystal is normally calculated by minimizing the free energy of the 
system subject to the condition that the crystal be electrically neutral 
throughout. However, as Lehovec (1953) has pointed out, the actual 
state of affairs is more complicated when we suppose that vacancies 
can only be formed at the surface of the crystal. (Compare also Frenkel 
1946, Grimley and Mott 1947 and Grimley 1950.) For definiteness we 
consider the case of sodium chloride. Here sodium-ion vacancies are 
more easily formed than are chlorine-ion vacancies. Hence when the 
temperature is raised from the absolute zero an excess of sodium-ion 
vacancies is emitted from the surface into the crystal, leaving a net 
positive charge on the surface. The resulting space-charge discourages 
the emission of further sodium vacancies and encourages tbe emission of 
chlorine vacancies. Calculation shows that in equilibrium the bulk of 
the crystal is electrically neutral, but that there is a positive charge 
on the surface, balanced by an equal and opposite negative charge cloud 
penetrating some distance into the crystal. Because of the presence of 
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this double layer there is a potential difference between the surface and 
interior of the crystal. 

It is generally admitted that vacancies can be formed at dislocations 
as well as at the surface of a crystal. Thus we should expect (Pratt 1957) 
dislocations in an ionic crystal to be charged (by having an excess of jogs 
of one sign), and to be surrounded by a sheath of vacancies of pre- 
dominantly the opposite sign. Clearly a grain boundary or mosaic 
boundary should also be charged and surrounded by a balancing layer 
of charge, whether we regard the boundary as essentially a junction 
between two separate blocks of crystal, or as an array of dislocations. 

We have pictured the charges on surfaces and dislocations as arising 
from the emission of an excess of vacancies of one sign. This is not 
necessary : any initial distribution of defects would re-arrange itself to 
give the same final equilibrium configuration. There should be similar 
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effects where Frenkel (interstitial plus vacancy) defects predominate, 
even though there is no difficulty in imagining their formation and 
disappearance in the body of the lattice as well as at surfaces and 
dislocations. 

Suppose that the temperature is low enough for the charge cloud around 
a dislocation to be immobile. The electrostatic field of the vacancies is 
such as to pull the charged dislocation back to the axis of the cloud if it is 
displaced, and so a certain minimum mechanical stress will be required 
to pull the dislocations away from their clouds and cause plastic flow. 
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Thus we have a mechanism which may in part determine the yield point 
of an ionic crystal. 

If a sodium chloride crystal contains divalent impurities mobile enough 
to play their part in forming the ion atmospheres the above account must 
be modified. At high temperatures there is qualitatively no modification, 
but at low temperatures the impurities cause the sign of the charge on a 
dislocation (or surface) to reverse. At one particular intermediate 
temperature the dislocations are uncharged, and at this temperature the 
yield strength of the crystal should drop to zero if the present mechanism 
were the only one determining it. Again, at low temperatures some of 
the impurity may precipitate out in a new phase, and thus reduce the 
concentration of impurity which take part in ion-cloud formation. If 
this is so the yield stress can have a succession of maxima and minima 
when plotted as a function of temperature. The measured points in 
fig. 1 do indeed show this sort of behaviour. 

The strength of ionic crystals will not, of course, be determined entirely 
by the mechanism discussed here. The factors commonly considered in 
connection with the strength of metals will also be operative, for example 
the elastic interactions between point defects and dislocations. We have 
thought it best, however, to present a preliminary account of the theory 
of charged dislocations ignoring these complications. Bassani and 
Thomson (1956) have dealt with the same topic, though from a different 
point of view. They showed that in NaCl a positive-ion vacancy is 
attracted into that part of the core of a (110) dislocation where the lattice is 
compressed, with a binding energy of about 0-4 ev relative to normal 
lattice positions. Divalent impurity ions were found to be little attracted 
into dislocation cores. From a knowledge of the binding energy of a 
point defect and a dislocation it is not possible to calculate a yield stress 
without introducing additional assumptions. Our approach, on the 
other hand, emphasizes the role of the charge cloud around a dislocation, 
and so we are able to obtain a yield stress by a purely electrostatic 
calculation. We may note that Bassani and Thomson’s model of a set 
of sites along a dislocation at which vacancies may be trapped with a 
certain binding energy is formally included in our ‘ vacancy source ’ 
model so long as the predicted charge on the dislocation does not require 
these sites to be completely emptied. Evidently the connection between 
the two models needs further examination. 


§ 2. THE CHARGE ON A CRYSTAL SURFACE 


We give first a brief sketch of the generalization of Lehovec’s results 
to the case where mobile divalent impurities are present. (A more 
detailed account will appear elsewhere.) 

Let N be the number of cation (or anion) sites and cNV the mean number 
of impurity atoms per unit volume. Let n,, n_, »; be the numbers per 
unit volume of cation vacancies, anion vacancies and impurity atoms at a 
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point where the electrostatic potential is v. (We suppose that v=0 at the 
surface of the crystal.) In what follows we assume that n,, n_, 2,, CN are 
all small compared with VN. We use an affixed co to denote the value of 
any quantity far from the surface of the crystal. 7’, e, k are the absolute 
temperature, the electronic charge and Boltzmann’s constant. 
If 7,., g_ are the free energies of formation of anion and cation vacancies 

the vacancy concentrations at a given point are 

n,=N exp {—(g,—en)/kT } 

n_=N exp {—(g_+ev)/kT}. 
The concentration of impurities must, according to Boltzmann’s baro- 
metric formula, depend on v through a factor exp (—ev/kT’). The 
constant corresponding to g,, however, is indeterminate, since for, the 
present we do not allow impurities to enter or leave the system. But at 
large distances from the surface n; must approach the bulk value cN, 
and. so 


n,=cN exp {—e(v—v,,)/kT }. 
Similarly n,* and n_© must agree with the values deduced neglecting 
the Lehovec effect. Thus 
n,2=N exp {—(g,—e0«)|kT}=Ne 
n_©=N exp {—(g_+ev,,)/kT} <> he eet a) 
=n -—n ~=— N(a—e) 
where 


a=g0-+ {($c)?-+exp [—(9,+-9_)/kT]}32. 2 tes Mecha 


A derivation of (2) has been given by Lidiard (1957). For our purpose 
it is enough to note that «=exp [—}(g,+g_)/kT] in the region of intrinsic 
conductivity but that at low enough temperatures the number of positive- 
ion vacancies equals the number of divalent impurities, so that z=c. 

vz, 18 determined by the first equation of (1) which may be put in the 
form 


ev./kT =n a+g,/kT. iy 5 bi 
Thus 
n,=Na exp p 
n_—=N(a—c) exp (—p) 
n,=Nc exp (—p) 
where 


p=e(v—v,,)/kT. 
The charge per unit volume, 


p=—e(n,—n_—n,)=2Na sinh p 


is related to the potential by Poisson’s equation V2v= —4zp/e where « is 
the static dielectric constant of the material, and so 


V2p= 2 sinh p 


eee was) 
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where 
= 8r7re*Na/ekT’. ori get demas dey Eee aa (2) 

The solution of (4) for the plane boundary has been given by Lehovec 
(1953). (It occurs also in the theory of colloids ; see for example Verwey 
and Overbeek 1948.) The potential has effectively risen from 0 to 
¥.» at a depth «~! below the surface. (Our « and Lehovec’s characteristic 
length A are related by «A= /2.) The charge Q on the surface and the 
field # just inside the crystal are given by 

ZexkT' nt Cv~ 6) 
= en op et a a 
The solution for an internal surface is clearly obtained by putting two 
surface solutions back-to-back. The charge on the internal surface is 


Pint=2Q. Sr ae ee MCT) 


§ 3. THe Force Required TO MOVE A CHARGED DISLOCATION 


We consider next the field about a dislocation. We should have to 
solve (4) in cylindrical coordinates. Unfortunately it is not possible to 
solve this non-linear equation exactly. If p is small enough for us to 
replace sinh p by p the solution is 

p=AK(xr)+Bl1y(«r) BP ap Ph emt (8) 
where Ky, J, are the modified Bessel functions normally so denoted 
(Watson 1945) ; the equation may also be solved if p is large enough for 
sinh p to be replaced by $ exp p. In any case it is not quite clear what 
boundary conditions we should impose at the dislocation line, and so we 
shall use the following argument to find the charge on a dislocation. 

Consider a mosaic boundary made up of edge dislocations spaced a 
distance h<x~! apart. Macroscopically the charge per unit area is 2Q. 
Microscopically, however, the charge is not uniformly spread, but resides 
on the dislocations. Thus the charge per unit length of dislocation is 


o=2Qh. ee le Oe de eet eet O) 


Let / increase. At first o increases linearly with h. But when / passes 
through the value «! the charge clouds of the dislocations will no longer 
overlap and the dislocations are effectively independent. Thus it seems 
reasonable to suppose that the charge on an isolated dislocation is given 
by (9) with h set equal to «1 times a numerical constant near to unity. 
If we take this constant to be unity the charge on an isolated dislocation is 


= VAN Pap peepee ah. (1) 

For lack of the precise solution of (3) we use (8) to tind the field surrounding 

the dislocation. If the field is to fall off with distance the term in /) must 
be rejected. The charge on the dislocation is then 

a dv <AkT 

og ST. diva) ai2e 
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using the relations 
K,;' @=—K,(2), lbmak (al 2. oS eee 


x>0 


Comparing with (8) we have 
V—V.,=(20/e)K (Kr). 

The potential due to the line-charge on the dislocation is 
—(2a/e) In r+-const. 

and so that due to the charge cloud surrounding it is 


Veloua=(20/e)[ Ky («r)+In («r)]+ const. 
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The restoring force on a charged dislocation. 


If the dislocation moves a distance 7 from the centre of the cloud, the 
cloud remaining fixed, it is subject to a force 


; d 
F=—o aS Ccloud 
= (20%/€)[A(1er) — (xr)3] 
by (11). If the movement is the result of the force bz due to the resolved 
shear stress 7 acting on a dislocation with Burgers vector 6 we must have 
F+6b7r=0 in equilibrium, or 
T=0-80(07«/eb) f(Kr) 
where f(%) (fig. 2) is the function x 1—K,(x) divided by its maximum 
value, 0:40, which occurs when x—1-02. The shear stress required to 
pull the dislocation away from its cloud is thus 


Y =7Tmax=0°8007«/eb. ° ° . ° ° (12) 
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We shall suppose that y, and g are represented adequately by the 
linear functions of temperature 


9.(T)=9,°—kT In A, 
g_(T)=g_°—kT In A_ 
so that in the intrinsic conductivity range 
a=(A,A_) exp {—4(9,°+9_/kT}. 


It is convenient to introduce in addition the quantity 


(13) 


a’—=exp (—g,/kT)=A, exp(—g,%kT). . . . (14) 


Evidently «’ is the concentration of positive-ion vacancies which we 
should infer if we neglected the presence of impurities and also the require- 
ment of bulk neutrality. We may call it the ‘ naive’ concentration. 

Figure 3 (a) illustrates how v,, varies with temperature. As a function 
of 1/kT' the curve In «’ is of uniform slope —g,°. The curve In « has the 
constant value In ¢ in the impurity conduction range, and changes over 
to a line of slope —3(g,°+g_°) in the intrinsic range. According to (3) 
ev,,/kT is numerically equal to the vertical distance between the In «’ 
and In « curves, and is positive where In«’ is uppermost. At a certain 
temperature 7’, the two curves intersect and v,, is zero. Following the 
usage in colloid theory we shall call such a temperature an isoelectric 
temperature. At 7’, the charge on a surface or dislocation (eqns. (6) 
and (10)) is zero and, according to (12), the yield stress Y vanishes as 
shown schematically in the lower part of fig. 3 (a). In fact at an iso- 
electric point the naive and true vacancy concentrations coincide and 
the neglect of electrostatic effects is justified. If 7’, lies in the impurity 
conduction range its value is 


Oe Ne). 


- Even in the absence of vacancies the two curves may intersect at high 
temperatures provided g,°<g_° and A,<A_. At such an ‘intrinsic’ 
isoelectric point g,=g_ (eqn. (13)), that is, the naive concentration of 
negative vacancies has just caught up with the naive concentration of 
positive vacancies. The position of the intrinsic isolectric point will be 
altered by the presence of impurities since « approaches 


bo-+{(de*+4,4_} 


rather than (A,A_)1/? at infinite temperature. 

We have supposed above that the impurity concentration c remains 
constant throughout the impurity range. In fact, if the heat of solution 
of the impurity is positive it will precipitate out in a new phase at low 
temperatures (provided equilibrium is attained) and we shall have 


c=A exp (—B/kT) Pe ee ou orl) 
with nearly constant A and B. As the temperature rises (15) will reach 


P.M. F 
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a value c=c, at which all the precipitate is re-dissolved and at higher 
temperatures c will remain at the value cy. Thus in the impurity range, 
at low temperatures, the In« curve will have a sharp knee and will 
intersect or fail to intersect the In «’ curve in one of the ways indicated 
in fig. 3 (b) according to the particular values of ¢,,, B, A,, A. Two of 
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these possibilities are shown in greater detail in fig. 3 (c), (d) together with 
qualitative sketches of the behaviour of the yield stress. In fig. 3 (d) 
there are two isoelectric points 7’, and 7’,, and a maximum in Y at 
T'2. (Since (12) is dominated by the behaviour of the hyperbolic sine 
the maxima of | ev,,/k7' | and Y will very nearly coincide.) In fig. 3 (c) 
there is a minimum in the yield stress at 7',,. (We ignore the dotted 


Charged Dislocations and the Strength of Ionic Crystals 83 
curves for the present.) The values 
T a=(B—g,)/k ln A 
TP as=9,/k In (cy?) 
T ga=T x= Blk In (Aley) 


are easily verified : the position of the maximum or minimum is inde- 
pendent ofg,. If there were a number of impurities with different values 
of A, B, cy the In« curve could in principle have a complicated zigzag 
shape and give rise to a succession of isoelectric points and maxima and 
minima. 

There is also the possibility that vacancies and impurities may associate 
at low temperatures (Lidiard 1954). As in the case of precipitation this 
will lead to a fall in the In« curve below a certain temperature, but 
the knee at 7’,, or 7',. will be smooth instead of abrupt, and the change 
of slope small. 


§ 4. THE VALUES OF g, AND g_ 


In most problems relating to Schottky defects in ionic crystals (con- 
ductivity, diffusion) it is sufficient to know only the sum g,+g_ of the 
free energies of formation, since it is this sum which determines the bulk 
concentration of defects. However, in order to predict yield stresses 
from eqn. (12) it is necessary to know the values of g, and g_ individually. 
Moreover, we have implicitly assumed that the values of g, and g_ for 
the formation of vacancies at a surface are still appropriate when the 
vacancies are formed at a dislocation, and this assumption evidently 
needs examination. 

It is clear that if dislocations are supposed to be in thermai equilibrium 
in a crystal g,+-g_ for a dislocation must be equal to g,-+-g_ for a surface, 
since otherwise Schottky defects would distil over from one to the other. 
We can establish the equality more directly with the help of an imaginary 
experiment. Consider a crystal containing, say, a single edge dislocation. 
We have to show that no work is required to transfer an Na—Cl pair from 
a jog on the dislocation to a kink in a step on the external surface. 
Cleave the crystal in two along a plane containing the dislocation line. 
This will require the expenditure of a certain amount of work. If the 
plane is suitably chosen one of the new surfaces formed will be flat, 
and the other will be flat except for a single step representing the termina- 
tion of the extra half-plane of Na—Cl pairs associated with the dislocation. 
Transfer an Na—Cl pair from a kink in this step (formerly a jog on the 
dislocation) to a kink in a step on the external surface. This requires. 
no work. If the halves of the crystal are now joined together again 
the work expended in cleaving is recovered. 

This argument cannot be applied to the transfer of a single ion, since 
the transfer is not a ‘repeatable step’ which leaves the physical state 
of the crystal unchanged. Thus we cannot infer that g, for the dislocation 


F2 
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is the same as g, for the surface. Strictly speaking even the transfer 
of a neutral pair shifts the position of a jog and alters the form of the 
dislocation and its relation to the free surface, so that the energies of 
cleaving and re-assembling do not exactly balance. The difference, 
negligible in the present context, provides a driving force which tends 
to make the dislocation climb out of the crystal. 

Etzel and Maurer’s (1950) experiments give the value 


a=5-2 exp {—1-01 ev/kT'} iy Me. eg 


for the concentration of Schottky defects in the intrinsic range. Thus 


(cf. eqn. (13)) 
9° +g_°=2:02 ev 


and A,A_=(5-2)?. It is reasonable to suppose that the frequency of 
vibration of an ion bordering a vacancy is reduced whether the vacancy 
is positive or negative. In this case both A, and A_ are greater than 


unity and 
Lage 2 oe, pte) Sg) eee ae 


This means that g, at room temperature can differ from g,=g_,° at 
0°K by a tenth of an electron-volt or so (compare (13)). 

The calculations of Mott and Littleton (1938) and their later modi- 
fications are mainly aimed at finding the sum g,°-+g_° which determines 
the bulk concentration of defects. The separate energies, g,’, g_’ say, 
required to remove a positive or a negative ion from the bulk of the 
crystal are found during the course of the calculation. When a positive— 
negative pair is replaced on the surface of the crystal or a dislocation 
an amount of energy L, the lattice energy per ion pair, is recovered, and 
60 g,°--g 99.49 '—L. 

In order to find g,° and g_® individually we have to decide what 
fraction of L is recovered when a positive ion alone or a negative ion 
alone is replaced on a surface kink or dislocation jog. At first glance 
we might think that this fraction would be the same for a positive ion 
as for a negative ion if only Coulomb forces and nearest-neighbour 
anion—cation repulsive forces were operative. This would give 
g,°=9, —3L, g.°=g_'—3L. This is Lehovec’s assumption. In this case 
it seems more reasonable to use the results of Mott and Littleton or of 
Brauer (1952), who neglect next-nearest neighbour interactions, rather 
than those of Bassani and Fumi (1954) who do not. Mott and Littleton’s 
figures give g,°—0:65 ev, g °=1-21 ev, Brauer’s g,°=0-70 ev, g_°=1-4]1 ev. 
However, the assumption that / is split equally in this way ignores the 
fact that on adding, say, a positive ion the ionic displacements in the 
vicinity of the kink or jog change from those appropriate to a negative 
kink or jog to those for a positive kink or jog : and conversely on addition 
of a negative ion. This means that the energy gained on adding a 
positive ion from infinity will be different from that gained on adding a 
negative ion, in the same way that the different displacements around 
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positive and negative vacancies make g,’ and g_’ different. Indeed it is 
possible that this difference is comparable with or even equal to the 
difference between g’ and g.’. In the extreme case g,° and g_° will 
each be equal to one half the energy of formation of a Schottky pair. 
Lehovec’s assumption corresponds to the opposite extreme where the 
displacements round a kink or jog are primarily determined by other 
factors and are not sensitive to the sign of the kink or jog. Evidently 
without detailed calculations all we can say is that probably g,° lies 
between 0-6 and 1-0 ev and that the values for surface and dislocation 
need not be identical. 
§ 5. Discussion 


As it stands the theory is unlikely to give good numerical results 
even with the correct values of g, and g_. The replacement of sinh p by 
p in (4) is not justified unless 

se a ee aon Bien LS) 


in the case of the plane solution. A similar or more severe restriction 
will apply to the cylindrical solution. With any reasonable values of the 
constants (18) is only obeyed within a few tens of degrees on either side 
of an isoelectric point. Even if the mathematical approximation were 
adequate outside these regions the high defect concentration which it 
predicts near a dislocation would require the physical assumptions to be 
refined. Moreover, in just those regions where the theory is at its best 
(near an isoelectric point) other hardening mechanisms, masked elsewhere, 
will prevent the yield stress falling to zero. However, we may hope that 
the experimental maxima and minima in the yield stress will coincide 
with the maxima and minima of | ev,,/k7' |. In the present paper we 
therefore limit ourselves to a tentative comparison of the critical 
temperatures in figs. 1 and 3. 

The points in fig. 1 show the yield stresses of a series of rock-salt 
specimens measured at various temperatures. The experimental details 
will be given elsewhere. Briefly, all specimens were cleaved from one 
large melt-grown single crystal, annealed for half an hour at 650°c, 
furnace-cooled to room temperature, and finally re-heated and held at 
the testing temperature for fifteen minutes before measuring the stress— 
strain curve. The divalent impurity content of the material is unknown, 
but it is estimated to be one part in 10°. For brevity we give the labels 
M,, M, to the maximum and minimum in fig. I. 

If g, is less than g_ than In «’ must be above In « at high temperatures. 
The transition from intrinsic to impurity controlled conduction must then 
be associated with a maximum in |ev,,/kT | (compare the dotted curves 
in fig. 3(c)). This could perhaps be identified with M, if allowance is 
made for experimental scatter and the uncertainty in the impurity 
content. Identification of the minimum M, with any of the possibilities 
of fig. 3 enables one to write 


94.=kT 4, In (1/c), 
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¢ being the free impurity concentration at M,. We then find that 7, is 
less than 0-6 ev unless (locating M, at 320°K) ¢ is less than about 10°%. 
Even if some of the impurity is precipitated or associated with vacancies 
such a value would seem unreasonably low. (Compare Haven’s (1955) 
measurements of the solubilities of Ca, Cd and Mn in NaCl.) This line 
of interpretation therefore leads us to a value of gy, considerably less than 
g_., and which is probably about 0-4ev. It is interesting to note the 
agreement with Bassani and Thomson’s (1956) value for the binding 
energy between a dislocation and a positive ion vacancy. 

An interpretation based on the alternative starting point that g, is 
greater than g_ cannot be sustained. For we should then have to suppose 
that the maximum M, was associated with the onset of precipitation 
and that the minimum M, represented a crossing of the In « and In ¢’ 
lines. For this to occur the heat of solution of the impurities would 
have to be greater than g,, i.e. greater than 1 ev, which would be quite 
inconsistent with the conductivity results. 

Thus we may say that comparison of our theory with experimental 
results on NaCl suggests that g, is considerably smaller than g_, perhaps 
only a quarter of it. It suggests that the maximum M, is associated 
with the change from impurity-controlled conduction to intrinsic conduc- 
tion, and that the minimum M, is associated with the onset of precipita- 
tion of the principal impurities. It would obviously be very interesting 
to have experimental results for crystals containing known amounts of 
deliberately added impurities. We should expect a reduction of strength 
at high temperatures (associated with the increased In «) with the removal 
of the maximum to higher temperatures (associated with a higher 
temperature for the ‘ knee ’ in the In « curve) and possibly the appearance 
of a hew minimum at intermediate temperatures. 

It is possible, however, that the fall in yield stress above M, has another 
explanation, namely that the charge cloud surrounding a dislocation is 
mobile enough to move with it at high temperatures. Electrical 
phenomena on deforming rock salt at various temperatures gives some 
support to this idea. 

Ifa dislocation moves away from its charge cloud an electric polarization 
is produced, proportional to the product of charge and distance moved. 
The movement of charged dislocations should thus give rise to a potential 
between opposite faces of a crystal, or to a flow of charge between 
electrodes connected by a low-resistance circuit (Pratt 1957). Plastic 
deformation is proportional to the distance moved by the dislocations, 
and so we might at first sight expect proportionality of voltage (or 
charge) and plastic deformation. However, though the dislocations are 
all of the same sign electrically, they can have opposite signs mechanically, 
in the sense that for a dislocation with a given Burgers vector we can 
always find another with an equal and opposite Burgers vector. Uniform 
deformation will involve the movement of nearly equal numbers of 
dislocations with opposite Burgers vectors in opposite directions, and 
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there will be no large net electrical polarization. In inhomogeneous 
deformation, on the other hand, a net flow of dislocations of one mechanical 
sign in the same direction can occur. Marked electrical effects of this 
_ type were indeed observed by Fischbach and Nowick (1955) under such 
conditions. 

Caffyn and Goodfellow (1955) observed smaller electrical effects in 
homogeneously deformed rock salt, and they found, moreover, that the 
effects decreased with temperature, disappearing altogether above 600°K. 
This suggests that at high temperatures the cloud follows the dislocation, 
so that no separation of charge is possible. In a temperature range 
where defects impeding a dislocation can only just keep up with its 
motion we may expect to observe a serrated stress-strain curve (compare 
Cottrell 1953). It is interesting to note that Classen—Nekludowa, (1929) 
observed jerky flow in rock salt at 500°K and above. 

We conclude by mentioning some phenomena which might be associated 
with charged dislocations. Though we have been considering Schottky 
defects there should be similar effects when Frenkel defects predominate. 

By plastic bending it is possible to introduce an excess of dislocations 
of one mechanical sign into a crystal. A crystal treated in this way 
should give large charge displacements or voltages when deformed homo- 
geneously. In such a crystal it would also in principle be possible to 
observe the converse effect, plastic deformation by the ‘ electrophoretic ’ 
movement of dislocations in an applied electric field. 

The movement of a dislocation under stress modifies the elastic constants 
of a material and, in an alternating stress field, there is a contribution to 
the mechanical damping. A charged dislocation should make analogous 
contributions to the dielectric constant and dielectric loss when it moves 
under the influence of an applied electric field. However, the restoring 
force which binds a dislocation to the centre of its charge cloud is pro- 
portional to r In (i/«7) rather than to r when it is given a small displace- 
ment r (compare fig. 2). Thus for a charged dislocation these effects, if 
observable, should be markedly non-linear. 

A dislocation is repelled by the charge on a crystal surface or internal 
mosaic boundary made up of dislocations with spacing less than «1. 
According to (6) and (10) the repulsive force near the wall is 7/0-80 or about 
four times the force necessary to tear a dislocation from its charge cloud. 
Thus when a dislocation has been pulled away from its charge cloud it 
may be held up by internal boundaries or the surface of the crystal 
until the stress is raised somewhat. This only applies when the disloca- 
tions in the boundary are closely spaced in comparison with x’. If 
«+ is small compared with the spacing each dislocation is surrounded 
by an individual ion cloud which completely screens its charge and 
electrostatically the boundary offers almost no hindrance to the passage 
of a dislocation through it. In a crystal with divalent impurities the 
screening length «~! drops sharply as we pass from the impurity to the 
intrinsic conductivity range. According to this model a mosaic boundary 
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behaves like a Venetian blind, closed to dislocations at low temperatures, 
but open at high temperatures. 

Again, two charged dislocations screened by their charge clouds exert a 
short-range repulsion on one another. If the dislocations are, say, a 
distance r=4x«-1! apart the electrostatic force between them is 80% of 
its unscreened value 202/er. Their elastic fields give rise to a force of 
the order of y.b?/27r, where p is a suitable elastic constant and 6 is the 
Burgers vector. These two forces are approximately equal if o is of the 
order of one electronic charge per atom plane. Consider, for example, a 
regular array of edge dislocations forming a vertical tilt boundary. If 
one dislocation, together with its ion cloud, is displaced out of the plane 
of the boundary the elastic fields of its neighbours tend to pull it back, 
but their electric fields tend to repel it further. In some circumstances 
an energetically reasonable compromise might be for successive disloca- 
tions to lie alternately to left and right of the mean plane of the boundary. 
Such an arrangement has been observed in ionic crystals (Mitchell 1956). 

Lehovec (1953) has shown that in a pure NaCl crystal there should be 
an extra conductivity proportional to the surface area, since in the ion 
cloud near the surface the concentration of cation vacancies (which 
earry most of the electrolytic current) is greater than it is in the body of 
the crystal. His analysis is easily extended to the case where there are 
mobile divalent impurities. Evidently the concentration of positive-ion 
vacancies near a surface exceeds or falls short of its bulk value according 
as V,, is negative or positive. Thus a surface increases the conductivity 
for temperatures where In «’ is above In « in fig. 3, but decreases it where 
In« is above Ina’. According to the argument leading to (7) and (9), 
one square centimetre of an internal boundary made up of dislocations 
spaced more closely than «~! should affect the conductivity in the same 
way as two square centimetres of external surface, while an isolated 
dislocation should be equivalent to a strip of external surface of width 2«~!. 

A number of crystallographic points need further discussion. Charge 
cannot form on pure screw dislocations. It may be energetically 
advantageous for a screw dislocation to turn towards edge orientation 
so that it may build up a charge with its balancing cloud. Again, if a 
charged segment of dislocation forming a Frank—Read source of length / 
expands into a large loop its original charge will not be spread uniformly 
over the loop. A simple argument suggests that only a segment of 
length / will be charged, and that the line joining its centre to the centre 
of the unexpended Frank-Read source will be parallel to the Burgers 
vector. Subsequent dislocations from this source will be initially 
uncharged. 
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ABSTRACT 


The decay of long-lived !°*Ho is discussed and results are given which 
suggest that the 817 and 718kev y-rays are each in competition with 
separate cascade pairs of energies 537 +-280 kev and 438 + 280 kev respectively. 
A level scheme is proposed consistent with the existence of a low lying 
rotational spectrum in the product nucleus 1°°Er. 


§ 1. INTRODUCTION 


TRRADIATION of holmium in the pile leads to the formation of two holmium 
activities attributed to the ground state and an isomeric state of 1*°Ho. 
The short-lived state of 27 hr half-life decays by 8 emission with the 
excitation of a number of states in 1°*Er including the first excited state 
at 80 kev: an exhaustive study of this mode of decay of Ho has 
been made by Graham, Wolfson and Clark (1955) and the state of 80 kev 
in 1°6Kr has been identified as 2+. The decay of the state of 1**Ho with 
a half-life greater than 30 years has been studied by Milton et al. (1955) 
who showed that the 6~ decay was accompanied by y-rays of energies 
845, 725, 282, 184 and 80 kev. Coincidence measurements suggested 
that these were emitted in a single cascade. 

The work described in the present paper was started as part of a 
programme of experiments on the alignment of the rare-earth nuclei, with 
a view to determining the magnetic moment of long-lived 1°*Ho and also 
to gaining further information on its decay scheme. In order to interpret 
such an experiment unambiguously, it is necessary first to know the 
decay scheme of the nucleus involved, including, if possible, the spins 
of the levels and the multipolarities of the y-rays (see for instance 
Blin-Stoyle and Grace 1957). A preliminary study of this isotope 
confirmed the presence of the five y-rays found by Milton et al. but 
led us to believe that the decay scheme differed significantly from that 
proposed by them. It was decided to use scintillation counters to 
investigate the y-ray spectra in coincidence and to extend the previous 
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work to include a study of summation peaks and of angular correlations 
so as to obtain an unambiguous decay scheme and if possible to limit 
the possible spin assignments of the levels. 


§ 2. Gamma-Ray SPECTRA 


A specimen of 200 mg of ‘ Specpure ’ Ho,O, was irradiated for 12 days 
in the A.E.R.E. pile at a neutron flux of 1-410! neutrons/cm?/sec. 
When the short-lived activity had been allowed to decay out a source 
of approximately one microcurie of 1*Ho was obtained : such a yield is 
concordant with that obtained by Milton et al. A spectrum of y-rays 
from this source in * good ’ geometry is shown in fig. 1. This was obtained 


Fig. 1 
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Spectrum of y-rays from 1°*Ho recorded with a 2in.x2in. Nal (TI) crystal. 
For meaning of symbols +, etc. see fig. 4. 


with the source 12 in. from a 2 in. diameter x 2 in. high Nal (TI) crystal 
(Harshaw) fitted to an E.M.I. type 6097 F photomultiplier. This instru- 
ment gives a resolution of about 8°% for y-rays of 840 kev. The spectrum 
was recorded on a 120-channel kicksorter of the Scarrott-Hutchinson 
type. The y-ray spectrum remained sensibly unaltered over the period 
-of a year and therefore if any contamination was present it must have 
been of half-life long compared with this. Although no chemical puri- 
fication of the source was carried out the close resemblance between our 
spectrum and that of Milton et al. (who did carry out a chemical separation) 
leads us to believe that our source was effectively pure. 

The energies of the y-ray peaks (y, etc.) are shown in the table and 
include the five y-rays already known. These energies were determined 
using standard sources of *4Mn, 1°7Cs and 23H» of energies assumed to 
be 840, 662 and 283 kev respectively. Also shown are the relative 
‘intensities which were estimated using calculated absorption coefficients 
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for Nal and graphical estimates of observed ‘ peak-to-total ’ values for 
these three sources. 

These results show clearly that the 280 kev y-ray is weaker than the 
other four. It should algo be pointed out that the weak peaks revealed 
at 438 and 537 kev cannot be attributed entirely to Compton effect : 
the broken curve shown under them represents an attempt to estimate 
the latter, using the data from the calibration sources. 

Because of the difficulty in interpreting the 280 kev y-ray as part 
of the main cascade, it was thought desirable to study the formation 
of ‘sum’ peaks to determine the total energy liberated in the decay 
chain. 


Gamma Rays from the Decay of 1®°Ho 


Our Milton et al.’s| Adopted Belative iene yaa 


Tahal value of value of value of y-Tay peaks 
“| y-ray energy| y-ray energy| y-ray energy Caled) Gang 
(kev) (kev) (kev) Ungated onan Mania 

WAL 80+43 80-3 +2 iy 80:3 +2 1-1 = — 
Vo 183-45 184+1 F 184-+1 1-0 1-0 1 
| 718 210) a Tob eES0 718 +10 1-0 OSes 
Sli Ao 845 +20 817-410 1-0 mae (Pie: 
ont OTT es 282 44 280-43 03+ | 0-154] 0-15¢ 
vom e4104815 = 438 (0.1L) |) (0:2), 
vy, | 528-415 =F 537 (0-1) Bein 


+ From f-spectrometer measurements by Graham and Clark (1955). 
t Not corrected for internal conversion of y;. 


The spectrum due to a 0-01 ue source in contact with a 4inx4in. 
crystal is shown in fig. 2. Sum peaks are seen corresponding to the 
main y-rays but there is no evidence for one corresponding to the simul- 
taneous capture of the 817, 718, 280 and 184 kev y-rays (the 80 kev 
y-ray contribution is expected to be difficult to observe because of its low 
energy and internal conversion). The intensity of this total energy peak 
is less than 2°%, of that predicted. Since both the ‘ gated’ spectrum 
described in §3 and the sum peak measurements were made with approxi- 
mately the same resolving time (2 psec) a finite lifetime to one of these. 
states could not account for the absence of a sum peak. We therefore 
conclude that the decay scheme cannot consist merely of a single cascade. 
chain. On the other hand it is clear that the 817, 718 and 184 kev 
y-rays are in cascade. 

Using the 4 in x 4 in. crystal in ‘ good ’ geometry, an attempt was made. 
to detect y-rays of energies greater than 817 kev. No peaks corresponding - 
to cross-over transitions were observed with an intensity greater than. 
1% of the 817 kev y-rays. 
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Fig. 2 
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Spectrum of y-rays from 1°*Ho showing sum peaks under conditions of poor 
geometry, recorded with a 4in.x4in. Nal (TI) crystal. 


Fig. 3 
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Spectrum of y-rays in coincidence with (A) ys, and (B) y,. 
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§ 3. CoINCIDENCE EXPERIMENTS 


The coincidence measurements were made with a pair of scintillation 
counters fitted with 2 in. crystals at a distance of 3 in. from the source. 
The output from one was fed through a single channel pulse height 
analyser, the output of which was used to gate a multichannel instru- 
ment which recorded the spectrum from the second counter. The 
results of spectra gated on the 718 and 817 kev peaks are shown in 
fig. 3: the intensities of the peaks are given in the table. Again it is 
evident that the 817 and 718 and 184 kev y-rays are in coincidence and 
that there is in each gated spectrum a y-ray of energy about 280 kev 
with about half the intensity of that in the ungated spectrum. Such 
a behaviour has been noted by Milton e¢ al. (private communication). 
In our gated spectra there is no significant difference between the energies 
of the two ‘ 280 kev’ y-rays observed. However, the fact that peaks at 
438 and 537 kev appear in coincidence with those at 817 and 718 kev 
respectively suggests that the 280+ 438 kev and 280+-537 kev groups 
are associated with the corresponding cross-over y-rays of 718 and 817 key. 

Angular correlation measurements showed that the 718-817 and the 
718-184 cascades were isotropic to within 1°, and 5°, respectively. For 
the 817-184 pair, measurements were made at 90°, 135°, 180°, 225° and 
270° and the angular correlation can be fitted by the function: 


W(0)=1+A,P.(cos 0)+A,P,(cos 6), 


with A,=—0-19+0-06 and A,=0-01--0-06. Such an angular correla- 
tion is not unique and does not permit an unambiguous assignment of 
spins to the levels involved. 


§ 4, Discussion 


The decay scheme of fig. 4 incorporates the idea of a pair of cascade 
transitions by-passing the 718 and 817 kev y-rays and appears to offer 
the simplest explanation to the observed discrepancies in the hypothesis 
of a straight cascade, namely : 


(i) The absence of a sum peak of a total energy of 2000 kev, although 
a strong peak of 1720 kev exists. 

(ii) The existence of y-rays at 438 and 537 kev whose total intensity 
approximates to that of the 280 kev y-ray. 

(iii) The apparent reduction in intensity of the 280 key y-ray in the 
“ungated ’ as compared with the ‘ gated ’ spectrum. 


The present evidence is insufficient to determine the sequence of 
the 718 and 817 kev y-rays, and similarly ambiguity exists as regards 
the intermediate level positions for the cascade transitions. Other 
possible explanations for these observations have been considered but 
the present scheme seems to offer the most straightforward key to the 
disintegration scheme. 
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If it is assumed that the transitions involve no higher multipoles than 
quadrupole then some limitations can be put on the level spins. From 
the angular correlation measurements of Milton et al., J, must be 4 and 
the energies of the y-rays strongly suggest that J, is the third member 
of the rotational sequence built up on this ground state with spin 6+. 
Since none of the angular correlations we have observed is of the type 
(J-+-4-— J+2-— J) we know that /, cannot be 6. In view of the transi- 
tions to the two members of the rotational sequence /,, 7, and not to 
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Proposed decay scheme for 1°*Ho (6-)!®Er, Energies shown are in kev. 


I, it seems probable that J, is 5. By admixing a small quadrupole 
constituent in the principally dipole transition, agreement can be obtained 
with either of the observed correlations with y,. Again, allowing a 
possible small quadruple admixture in ys, the observed correlation 
Yo-y4 is not consistent with the unique result which must occur were 
I, greater than 6. Approximate equality in the excitation of the different 
states I,, I,, I, shows that the majority of the transitions proceed 
through the state 7,. In order to account for this, it is suggested that 
both 16*Ho and the level J, possess odd parity and that the parity of the 
lower states in !*Er are all even. 
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ABSTRACT 


Helium atoms have been injected into copper by alpha-particle bombard- 
ment in a cyclotron. On subsequent heating, the helium precipitates out 
as small gas bubbles in those parts of the metal where the helium atoms 
receive vacancies from nearby sources. From the distribution of the gas 
bubbles it is concluded that grain boundaries are effective sources of 
vacancies, but twin boundaries are not. Isolated sources of vacancies 
exist within the crystals, but there is no evidence that the dislocation net- 
work generates appreciable numbers of vacancies. 


§ 1. INTRODUCTION 


WHEN the temperature of a crystalline solid is increased the equilibrium 
number of vacancies in it increases steeply. In most metals these 
vacancies cannot originate by vacancy-interstitial pair creation in the 
perfect lattice because the energy of formation of an interstitial is much 
higher than that of a vacancy. The vacancies must enter good regions 
of the crystal from various sources, for example, grain boundaries, sub- 
grain boundaries, dislocations, and free surfaces, where it is not necessary 
to form interstitials with vacancies. Evidence for such sources has 
previously been indirect, depending for example on observations of voids 
formed during the interdiffusion of two metals (Balluffi and Alexander 
1952, Barnes 1952, Biickle and Blin 1952, Seith and Kottmann 1952). 
In the experiment described below we have developed a method for 
locating vacancy sources in a metal and for observing, on a microscopic 
but not atomic scale, these sources in operation. Helium atoms are 
first injected into the metal by bombardment with alpha-particles. On 
subsequent heating, these atoms attempt to precipitate within the metal 
in the form of gas bubbles, and to acquire the extra space necessary for 
this they capture large numbers of vacancies. Thus a blanket of bubbles 
forms in the regions of the metal immediately accessible to vacancies 
from the sources, but until these bubbles have become sated with 
vacancies, no bubbles are formed in the remaining regions, which they 


blanket. 
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§ 2, MerHop AND RESULTS 


Johnson—Matthey spectroscopically standardized copper was bom- 
barded with 30 Mev (approximately) alpha-particles from the Birmingham 
cyclotron until 2-5 kw hr had been dissipated in it. This is estimated 
to have deposited helium to a concentration of 1 atm%, in a band of 
metal about 0-003 em thick lying at a depth of 0-12 cm below the 
bombarded face of the specimen. Being brazed to a water-cooled 
copper target holder, the specimen did not heat above 200°c during the 
irradiation. After allowing the radioactivity to decay, pieces were 
cut from the specimen and individually heat-treated. Each piece was 
then polished, mechanically and electrolytically, and examined micro- 
scopically. 

Those pieces not heated after bombardment showed no bubbles. Only 
after heating to 650°c for 1 hour did bubbles first appear, and then only 
in localized regions. Figure 1, Pl. 4 shows a section cut perpendicularly 
through the band of helium, taken from a specimen heated at 750°c for 
1 hour in vacuo. The gas bubbles are localized in regions which appear 
dark and which occur (a) where grain boundaries cross the band, (b) on 
the periphery of the band, and (c) as small circular islands within the 
band. Electron micrographs of two-stage carbon replicas taken from 
the surface show that individual bubbles are less than 2 x 10~° cm across. 
Figure 2, Pl. 4 gives an example showing several islands of bubbles. It 
will be seen that the bubbles appear to be crystallographic in shape, but it 
is not yet clear whether this is their natural shape or is a consequence of 
the electropolishing treatment. 

Prolonged heating, particularly at high temperatures, causes the 
bubbled zones to spread farther out from the vacancy sources until 
eventually the whole band is uniformly filled with bubbles. This 
uniformity of the final distribution shows that the process does not 
occur by the migration of helium atoms to the sources of vacancies. 
The rate of advance of the bubbled zones is greatest near the periphery 
of the band, where the concentration of helium is lowest, and least in 
the centre of the band, where most of the helium is deposited. This is 
to be expected since, for a given rate of generation of vacancies from a 
source, the width of the bubbled zone fed by that source should increase 
at a rate inversely proportional to the concentration of helium present. 

These various effects have been observed repeatedly in copper speci- 
mens. Preliminary experiments have shown closely similar effects to 
exist also in beryllium. 

Figure 3, Pl. 5, taken from a copper specimen heated to 700°c for 2 hours, 
is a section almost parallel to the plane of the band. Such a taper section 
magnifies the width of the band and is more informative. Many coherent 
twin boundaries are visible and it will be seen that they do not assist the 
precipitation of the helium atoms. Non-coherent twin boundaries are also 
visible, particularly in the central grain of the photograph, and it will be 
ceen that these also, when not connected to grain boundaries, do not assist 
the creation of bubbles, Some single-ended sources of vacancies attached 


Observation of Vacancy Sources in Metals 99 


to grain boundaries are also visible. It is thought from these, and from 
similar examples seen in other specimens, that they are non-coherent 
twin boundaries attached to grain boundaries; and thus that these 
non-coherent boundaries can accelerate the transport of vacancies but 
not create them, whereas coherent boundares do not accelerate the 
transport of vacancies. 

§ 3. Discusston 

The conclusion that the bubbles grow by the capture of vacancies is 
based on the peculiar distribution of the bubbled regions in incompletely 
annealed specimens. Particularly striking is the observation that 
bubbles form along the edge of the helium band before they form in the 
centre, despite the concentration of helium being highest in the centre. 
There is no hindrance to the formation of bubbles along the edge, since 
these outlying regions of the band are accessible to vacancies arriving 
from all sources outside the band. 

Only two major sources of vacancies within the metal have been 
positively identified in these experiments ; grain boundaries, and certain 
sources within the grains which, because they produce bubbled zones 
that are practically circular in all sections taken through the band, 
appear to be point sources. The nature of the latter is not clear. 
They may be small cavities, present in the specimen before heating, 
which lose vacancies and close up during heating. Or they may be 
points, such as dislocation nodes, that can emit vacancies provided these 
vacancies can reach them along suitable channels connecting them to 
grain boundaries. Certainly, dislocation networks in the crystal, if 
acting at all as vacancy sources, are not doing so to a significant extent 
along most of their length. 
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ABSTRACT 


This note gives a brief mathematical treatment of the relationship between 
the transient stress in a semi-infinite visco-elastic rod and the attenuation 
coefficient and phase velocity as functions of steady state frequency. 
Except in the case of a delta function stress input, the expression obtained 
disagrees with Kolsky’s result (1956). The reason for the failure of the 
latter is pointed out. 


Unper certain conditions, which are usually satisfied, the general solution 
for the passage of a stress pulse along a semi-infinite visco-elastic rod can 
be related to the steady state solution. We use the result obtained by 
Berry and Hunter (1956) by means of Laplace transform methods that 
the stress at time ¢ at a point distant x from the origin is 
I) (pare 
Gf, t)=— xa a(0, s) exp [—p(s)a-+st] ds, x>0, t>0, (1) 
y—tco 
where ~ 20 
FLGo= | o(0, t) exp (—st) dt 
0 


and y.(s) is determined from the equation 


pps? | g(t) exp (—8t) di =a) a ee ee 


so that it is positive when s is real and positive. In this expression p is 
the density of the rod and g(t) is the strain per unit stress at time ¢ after 
the application of a constant stress. We suppose that the rate of creep 
at constant stress tends to zero as f—> oo so that g(t) can be represented 


by the expression 2 
=a | G@jil-exp (stride: a cn 


where G(7) is an unnormalized distribution function of retardation times 
of creep (Gross 1953). Substituting this form of g(t) into eqn. (2) we find 


that G(r) der 
piA== ps" Ga | a] pee ee ee ie) 
and it follows that p(s) is regular in the half-plane R(s)>0. 
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We limit the class of stress pulses applied at the origin to 5-function 
pulses and finite pulses of finite duration, ¢,, say. In the first case o(0, s) 
is a constant, while in the second it is of the form 


ai [°F exp (—at at ee eas) 

and is thus regular in ®(s)>0. We may then put 
o(0,4a)=A(a)—tBlw)so we (BY 
where Alois {10 coswtdt, B(w)= I, f(t)sinwtdt  . . (7) 


A(w) and B(w) being even and odd functions respectively (constant and 
zero for a 5-function). 

Under the above restrictions the integrand in eqn. (1) is regular in 
X(s)>0 and y may be taken as zero. It can be shown from eqn. (4) that 
the functions «(w), c(w) defined by 


a(w)=R[w(tw)],  c(w)=e/B[mliw)} . - - - + (8) 
are both even, positive functions of w. Then, introducing the functions 


(7) and (8) and making a change of variable, the integral (1) may be 
reduced to the real integral 


o(2, j= - ii [A(w) cos w(t—ax/c)+ B(w) sin w(t—az/c)|exp(—ax)dw. (9) 
at A) 


To show that «(w) and c(w) are the attenuation coefficient and phase 
velocity, respectively, we consider the stress in the rod due to the im- 
position of the sinusoidal stress 

a(0, t) =a, exp (dwt) 
at the origin. Then 
a(0, s)=a,/(s—iw) 
and, from eqn. (1), we have the result 
spo 7 ee eT a, een (10) 


ae O= Q71 S—tw 


y—to 
When ,(s) is defined by egn. (4), it is not difficult to show, by means of 
the residue theorem, that as t->0o the stress given by eqn. (10) tends to 
the steady state value 
o,(x, t) =o, exp (—ax) exp [tw(t—2/c)]. ay Oe ee Le 

Kolsky (1956, eqn. (5)) has given an expression for a stress pulse in a 
visco-elastic rod which consists of twice the first term in the integral (9). 
Representation of the stress input by a cosine integral alone is not valid 
since it is then, in effect, an even function of ¢ and a second pulse is implied, 
originating in a negative interval of time. In eqn. (9), however, half of | 
the applied stress is represented by a sine integral which, being an odd 
function of t, makes the sum of the two terms zero for <0. When there 
is a 5-function stress at the origin the two expressions agree so long as the 
§-function is defined symmetrically when it is subjected to the Fourier 
cosine transform and unsymmetrically for the Laplace transform (Sneddon 
1955). 
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It should be noted that the phase velocity c(w) is a function of both the 
real and imaginary parts of the complex modulus H*(iw) (Gross 1953). 


It is a consequence of eqn. (2) that 
2 


aie ee lb tl” 
E*(iw)=Hy(w)+1h (w)= [w(iw) |? ? ( ) 
and it follows by means of eqn. (8) that 
qa IEEE). S) oe 
pH, + /(#Y+#,’)] 
The velocity of propagation of the wave-front, c,, is readily obtained by 
considering the integral (1), which is zero if 
R(st—px) > — co as |s|— oo in | arg s|<$z. 
We find, using eqn. (4), that 
seeps B (1/8) = v/ (pa) 


or, from eqns. (8) and (12) 


EB , 
ee(co)=, / AX) eres 
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CORRESPONDENCE 


Propagation of Dislocations in Electron Irradiated 
Lithium Fluoride 


By A. D. WHAPHAM 
Metallurgy Division, A.E.R.E., Harwell, Berks. 


[Received October 14, 1957] 


A stupy of the properties of dislocations in irradiated materials is 
fundamental to an understanding of the effect of irradiation on mech- 
anical properties. The development in recent years of etching techniques 
which produce etch pits sited at the point of emergence of dislocation 
lines at the surface of a crystal has provided a useful method of observing 
dislocations. The 50% hydrofluoric acid, 50% acetic acid and ferric 
fluoride etch for lithium fluoride, developed by Gilman and Johnston 
(1956) and used as a tool in the present investigation, is particularly 
good since it reveals all dislocations, does not seem to depend on impurity 
atoms and does not require any annealing of the specimen. 

Dislocations have been produced in a Harshaw lithium fluoride (100) 
cleavage face by micro-hardness indentations both before and after 
irradiation by 1 Mev electrons at a flux of 34acm~? for 1 hour. 
Differences in dislocation behaviour were revealed by subsequent etching 
and the appearance of the cleavage face after treatment is shown in 
the figure (Pl. 6). The etching technique worked reliably on the irradiated 
material. Identations 1 to 4, made before irradiation with a load of 
about 2g on the Vickers diamond indenter, produced star patterns of 
dislocations with arms radiating along the ¢100) (screw dislocations) 
and <110) axes (edge dislocations). (These star patterns are essentially 
the same as those etched before irradiation (not illustrated) showing 
that no appreciable movement of these dislocations has occurred during 
the irradiation.) The edge dislocations have been propagated by the 
stress field to about 4:5 diameters and the screw dislocations to about 
3 diameters of the central impression. In the case of indentations 
5 to 10 made after irradiation (5, 6, and 7 with 2 g load and 8, 9 and 10 
with 10g) the edge dislocations have only been propagated to about 
1 diameter from the centre of the impression while the motion of screw 
dislocations has been so restricted that no evidence of propagation is 
seen on the photograph. 

It is known that irradiation produces hardening in most materials. In 
the case of alkali halide crystals this has been shown by Podachewsky 
(1935) for irradiation with ultra-violet light and x-rays and confirmed 
in the case of x-irradiation by Pratt (1951). It has been shown also 
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for fast proton irradiation (Vaughan et al. 1953) and for x-irradiation 
and 1 Mey electron irradiation (Westervelt 1953). In metals neutron 
irradiation increases the tensile strength and hardness (cf. McReynolds 
et al. (1954) in the case of aluminium, Gieb and Grace (1952) for 
molybdenum, Kunz and Holden (1954) iron, zinc, Meyer (1954) mild steel 
and Billington and Siegel (1950) for stainless steels). This hardening can 
be due either to difficulty of generation or of propagation of dislocations. 

In the present work indentation has generated dislocations in the 
irradiated material (large numbers of dislocations closely clustered around 
identations 8, 9 and 10 can be clearly seen), but no information about the 
generation stress is given by the experiment. The experiment proves, 
however, that the stress for propagation of the dislocations is raised by 
irradiation since the dislocations produced in irradiated crystal do not 
glide away from the indentation as readily as those produced before 
irradiation. Hence irradiation produces lattice friction. The origin of 
this friction is not known but it may be due to isolated defects or aggre- 
gates distributed in the crystal lattice. 

It is understood that Gilman (1956) has made a similar study of the 
propagation of dislocations in neutron irradiated lithium fluoride but, as 
pointed out by Lomer (1956), complications arise in the interpretation of 
the results since nuclear disintegration of the lithium occurs when it is 
bombarded with neutrons. (No such reaction can occur with 1 Mev 
electrons.) 

At present, work is in progress to determine the magnitude of the 
stress necessary to generate and propagate dislocations in lithium fluoride 
subjected to various amounts of electron irradiation. 
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The Half-Life of Thallium 196 and the Multipolarity of the 
426 keV Transition in Mercury 196 


By G. ANDERSSON, J. O. BuR@MAN and B. JuNG 
The Gustaf Werner Institute for Nuclear Chemistry, Uppsala, Sweden 


[Received October 4, 1957] 


Ty a b-spectrometric investigation of neutron deficient Pb and Tl isotopes, 
Andersson et al. (1955) observed the conversion lines of a 426 kev y-ray. 
This was identified as a transition in '°°Hg, earlier reported by Staehelin 
(1952), who studied the decay of !°*Au. Further arguments for the 
assignment were obtained from excitation data and from the non- 
appearance of the rather strong K-line in the conversion electron spectrum 
of an electromagnetically separated mass 197 sample. From a number 
of successive automatically recorded spectra the half-life, which should 
thus be characteristic of 1°°Tl, was very preliminary estimated as ~4 hr. 

Later the growth and decay of the K-line was more accurately measured 
(Andersson ef al. 1957). From a straight line over several half-lives 
the longer-lived component was determined as (2:4+-0-2) hr. Its growth 
indicated a 35 min parent, which was also consistent with the mass 
assignment 196, as 1°®Pb had independently been assigned a 37 min period. 
Knight and Baker (1955), who used enriched Hg targets to produce 
Tl-activities, report 3 hr half-life for 1°°Tl in acceptable agreement 
with the value quoted above. 

According to Staehelin (1952) the 426 kev transition de-excites the 
lowest excitation level of 1°*Hg. It should thus be of pure E2 character, 
which was supported by the measured conversion coefficients. The 
K/L conversion ratio obtained by Andersson ¢¢ al. (1957), however, 
indicated a considerable M1 admixture. 

We have now reinvestigated the pertinent conversion lines at about 
0-3°,, resolution in a double-focusing beta spectrometer (Arbman and 
Svartholm 1956) using mass separated sources (cf. Andersson 1957). 
The half-life of the K-line in 4d=196 was measured as (1-8-+-0-1) hr and 
the K/L ratio as 2-75-- 0-28 in good agreement with the value 2-86 obtained 
for electric quadrupole radiation by interpolation from Rose’s tables 
(1954-1956). 

The cause of the earlier discrepancy was found to be an M1 y -ray in 
the decay of 2-7 hr 19’T1, coinciding energetically with the '°*Hg transition 
within the resolving power of the beta spectrometer. In the above- 
mentioned measurement on A197 (Andersson et al. 1955) the sample 
must have been just weak enough for the K-line not to appear. It should 
be especially pointed out that 197T] grows from 42 min 1°’™Pb, that is, 
with about the same parent half-life as !°°Tl, which explains the result 
reported by Andersson et al. (1957). 


P.M. : 
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The present investigation illustrates the difficulties encountered in 
nuclear spectroscopic work on mixtures of several nuclides and emphasizes 
the importance of mass separated samples. 
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REVIEWS OF BOOKS 


Light Scattering by Small Particles. By H.C. vax pe Hutst. (London: John 

Wiley and Sons. Inc.) [Pp. xiii+-470.] 96s. 

MANY current scientific investigations depend upon the observation and inter- 
pretation of the scattering of electromagnetic radiation by small particles. 
These range from particles of colloidal metal in solid and liquid suspension 
through spherical water particles and ice crystals in the atmosphere to the 
particles of planetary atmospheres and of the dust clouds of interstellar space. 
In all these investigations the Mie theory based on the rigorous solution of 
Maxwell’s equations for the scattering of a plane wave by a homogeneous sphere 
is of primary importance. There is, however, a wide gulf between the formal 
solutions of the equations, which lead to complicated formulae, and practical 
applications ; the average experimental worker has difficulty even in applying 
the mass of tabulated numerical data which is now available on account of the 
pitfalls of interpolation. 

Professor Hulst gives an exhaustive treatment of the Mie theory for scattering 
by spheres over the whole range of sizes. The functions which are most needed 
for the interpretation of experimental observations on Raleigh scattering, 
Rayleigh—Gans scattering, scattering by large dielectric spheres and by reflecting 
metallic spheres are presented in graphical and in tabular numerical form. A 
detailed account is given of the best methods to be followed for obtaining 
numerical results and valuable suggestions are made for reliable interpolation 
in the published tables. Throughout the book emphasis is laid on the dis- 
cussion of the physical significance of the mathematical results. 

Although the book is mainly concerned with the Mie theory, there is an 
adequate account of other treatments and of light scattering by particles of 
non-spherical form. An interesting review of the author’s own work on 
atmospheric phenomena including the theory of the glory or Brocken spectre is 
included. 

Part II of the book, which will be of great value to physicists, meteorologists, 
astronomers and chemists, is concerned with practical applications. There is a 
useful discussion of the properties of particles which may be investigated by 
observations of light scattering, of the measurements which should be made, 
and of the information about the scattering particles which can be deduced from 
them. J. W.M. 


High Energy Physics : Proceedings of the Seventh Annual Rochester Conference. 
Edited by R. G. Sacus. (New York: Interscience Publishers Inc.) 
[Pp. 482.] $4.50 

THe annual Rochester report is an invaluable guide to the present state of high- 

energy physics in the days of a literature too vast for anyone to cover. The 

introductory surveys to the eleven sessions take up over a third of it, and 
enable any nuclear physicist to understand the current problems and the work 
towards their solution presented by the specialists in the remainder of the 
report. This year’s report comes when there is a pause in the flow of the subject. 
as last year’s did not, and is thus more useful as a reference work. Electro- 

_ magnetic investigations into the structure of the nucleon, the role of the dis- 

persion relations in pion physics and the extent of their mathematical vigour, 

the phase-shift analysis of nucleon-nucleon scattering and our understanding 
of the parity puzzles can now all be summarized in a reasonably compact form. 

Data are also being steadily accumulated in the less familiar fields of extremely 

high-energy cosmic rays and reactions involving antiprotons and strange 

particles, particularly on the interaction between nucleons and K mesons. 
A change of typewriter has improved the legibility, and photographs and more 

complicated diagrams can now be reproduced. D.C: 
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The Detection and Measurement of Infra-red Radiation. By R. A. SMITH, 
F. E. Jones, and R. P. Crasmar. (Oxford University Press, 1957). 
[Pp. xiii+458.] 70s. 

Tue major part of this book (about 250 pp.) is concerned with the detectors of 

infra-red radiation. Particular care has been taken to relate the actual per- 

formance of the various types of thermal and quantum detectors to ultimate 
limits set by fundamental physical principles. For the first time the science 

(as opposed to the art) of infra-red detection has been presented in book form 

and this is likely to remain the standard work on the subject for a considerable 

time. 

Another main section (about 140 pp.) deals with other experimental 
techniques used in the infra-red region. Here the treatment is uneven. Some 
chapters are very good but, for example, that on the design and performance 
of infra-red spectrometers is rather brief. It is a pity that no theoretical 
account has been attempted of the performance of a spectrometer as a@ whole in 
the case—common in the infra-red region—of an energy limited instrument. 
The reader of this book could not then have been left with the impression 
(possibly unintended) that diffraction grating spectrometers are used only when 
the higher resolution is required. Many are now being produced commercially 
which will mostly be used to give higher intensity spectra (better signal—noise 
ratio) than the more usual prism instruments. 

The publication of this book seems to have been rather slow for important 
and not very recent topics such as the indium antimonide photoconductive cell 
and F-centre filters are not mentioned. However, these are minor points of 
criticism. This book is an important addition to the infra-red literature. 
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BOOK NOTICES 


Solid State Physics. Advances in Research and Applications. Vol. 4. Edited 
by Freperick Serra and Davin TurNnButy. (London: Academic Books 
Ltd.) [Pp. xiv+540.] $12.00 
The fourth volume of this useful series contains the following articles : 


“ Ferroelectric and Antiferroelectrics > by Werner Kinzig (General Electric 
Research Laboratory, New York). 

“Theory of Mobility of Electrons in Solids” by Frank J. Blatt (Michigan 
State University). 
“The Orthogonalized Plane-Wave Method” by Truman O. Woodruff 
(General Electric Research Laboratory, New York). 

“ Bibliography of Atomic Wave Functions ” by Robert 8. Knox (University 
of Rochester) ; and 

“Techniques of Zone Melting and Crystal Growing” by W. G. Pfann (Bell 
‘Telephone Laboratories). N. F. M. 
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(b) < 1000 


Etched zine single crystal—surface parallel to basal plane. 
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Fig. 1 


X-ray diffraction pattern of mildly rolled polyethylene film. X-ray beam 
parallel to direction of rolling. Plane of rolling horizontal. 


X-ray diffraction pattern of strongly rolled polyethylene film. X-ray beam 
parallel to direction of rolling. Plane of rolling horizontal. 
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Fig. 3 


X-ray diffraction pattern of strongly rolled and subsequently lightly annealed 
polyethylene film. X-ray beam parallel to direction of rolling. Plane 
of rolling horizontal. 


Fig. 4 


X-ray diffraction pattern given by a polyethylene strip which was first. fully 
drawn, then annealed at 110°c and subsequently redrawn ~60%. —_Direc- 
tion of drawing vertical. 


R. S. BARNES et al. Phil. Mag. Ser. 8, Vol. 3, Pl. 4. 


Section perpendicular to the helium band in copper after an anneal of 1 hour 
at 750°c. 500. 
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Electron micrograph of a two-stage carbon replica taken from the surface of 
copper containing islands of helium bubbles. x 7000. 
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Taper section of the helium band in copper after an anneal of 2 hours at 700°c 
(same surface as seen in fig. 2). 500. 
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Etched (100) cleavage face of the electron irradiated lithium fluoride specimen. 
The cube axes are horizontal and vertical (x 500). 


